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SUMMARY 


Our experimental and theoretical program has been focused on the 
growth and characterization of GaAs. In our crystal growth work we have 
pursued the study of quantitative relationships between 'growth parameters 
and properties. In our characterization studies we have investigated the 
determination of the compositional, structural and electronic properties of 
GaAs on a macro- and microscale. We have developed a. detailed theoretical 
model of electroepitaxial growth and we have found it to be in excellent 
agreement with experimental results. The growth rate and composition of epi- 
taxial layers can be controlled by current density. Electroepitaxial layers 
are structurally and electronically superior to those grown by thermal liquid 
phase epitaxy. The growth rates attained by electroepitaxy thus far are one 
to two orders of magnitude greater than those obtained by standard liquid 
phase epitaxy and only one order of magnitude smaller than those obtained in 
melt growth. In addition to electroepitaxy we have completed the design, 
construction and testing of apparatus for the growth of GaAs crystals from 
the melt. Regarding characterization, we have developed novel techniques for 
obtaining two-dimensional microprofiles of- the carrier concentration in GaAs 
(sc ann ing IR absorption) and microprofiles of the minority carrier characteris- 
tics (SEM-EBIC mode). In addition, we have .developed derivative surface photo- 
voltage spectroscopy, a new' method, which permits the simultaneous determina- 
tion of critical point transitions and transitions involving levels in the energy 
gap. . At the same time we are nearing construction of standard characterization 
techniques such as transient deep level spectroscopy, photoluminescence and 
cathodoluminescence spectroscopy. Our aim is to utilize our novel and selected 
standard characterization techniques for investigating the electronic properties 
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of GaAs,. such as carrier concentration, compensation ratios, minority carrier 
lifetime, and diffusion length on a microscale, and for relating these proper- 
ties to materials parameters (compositional and structural) which in turn we 
will relate to growth parameters. Such relationships should permit the optimiza- 
tion of GaAs growth under zero-gravity conditions as well as the assessment and 
exploitation of the potential of GaAs in device applications. 

INTRODUCTION 

It is now generally accepted within the scientific and engineering community 
that to insure successful processing of materials (e.-g. , crystal growth) in 
space extensive work must be carried out on the ground. Ground-based work 
must be aimed at establishing the experimental parameters and conditions for 
the optimization of the unique advantages of zero gravity conditions. In addi- 
tion, experimental techniques and theoretical understanding must be developed 
which, permit the reliable characterization of the material (before and after 
space processing) and, thus, a detailed analysis of the implications of zero 
gravity environment to the advancement of basic science and’ to the realization 
of unique technological applications. Actually, the results obtained from the 
growth of indium antimonide in the Skylab mission and from the growth of ger- 
manium in the ASTP mission surpassed all expectations only because both systems 
were investigated extensively for many years prior to the design performance and 
analysis of those experiments. 

Accordingly, our ground-based research effort is aimed specifically at 
the optimization of space processing of GaAs by (a) developing an experimental 
and theoretical framework for the detailed compositional, structural and 
electronic characterization of GaAs on a macro- and microscale, (b) developing 
a quantitative understanding of the relationships between growth parameters and 
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the properties of GaAs, and (c) assessing the key GaAs parameters limiting 
device performance. 

In this second annual report on "Crystal Growth of Device Quality GaAs 
in Space" we present a summary of our results on the crystal growth and on 
the characterization of GaAs. In addition, we present more detailed discus- 
sions of some of our results as they have recently appeared or will soon 
appear in the open literature. 

Under "Proposed Work" we. discuss our proposed research plans for a one- 
year period beginning April 1, 1979. Some new research areas on growth and 
characterization will be experimentally initiated, whereas work now in progress 
in solution growth and in characterization will be further pursued. 

CRYSTAL GROWTH 
Solution Growth 

Thermodynamics predicts the solution growth., requiring lower growth 
temperature than melt growth, should yield single crystals of far better 
quality than those obtained from melt growth. Actually, in the case of GaAs, 
using liquid phase epitaxy (LPE) the quality of the. layers and their suitability 
for device fabrication and device performance is essentially limited by their 
interactions with the substrate materal which is generally of poor quality. 

The standard LPE technique, because of its very low growth rates, cannot be 
realistically considered for the growth of bulk crystals which could replace 
those grown from the melt. 

Electroepitaxy does not in principle have the growth rate limitations of 
standard LPE. In electroepitaxy crystal growth takes place by passing electric 
current through the growth interface while the temperature of the overall 



system is maintained constant. The required supersaturation for crystal growth 
is brought about by electromigration of solution species to the growth inter- 
face under the applied electric field. Thus, the slow diffusion-controlled 
transport of solute which controls the growth rate in standard LPE is of no 
consequence in electroepitaxy where significantly faster transport of solute 
takes place through electromigration in an electric field. 

We have completed the first stage of extensive experimental and theoretical 
studies which have resulted in the development of a comprehensive model of 
electroepitaxy. Thus, on the basis of mass transport principles a theoretical 
model of electroepitaxial growth — current-controlled LPE — was developed which 
defines explicitly the contribution of the Peltier effect (at the solid-solution 
interface) and that of solute electromigration to the overall growth process. 
According to the model the contribution of electromigration to growth is dominant 
in the absence of convection in the solution, whereas the contribution of the 
Peltier effect can be dominant in the presence of convection. On the basis of 
the model expressions were derived which relate quantitatively the growth velocity 
to growth parameters. The model was found to be in excellent agreement with 
extensive experimental data on the electroepitaxial growth of GaAs from Ga-As 
solution. The detailed treatment of this model and of related experimental 
findings is presented in the attached recent publication on "Liquid Phase Epitaxy: 
Growth Kinetics". 

In this report we will present typical recent results illustrating the 
uniqueness and the potential of this technique. We have found that electroepitaxy 
has significant advantages over' thermal epitaxy. (1) The layers grown by 
electroepitaxy exhibit a far superior surface morphology (Fig. la) than those 
grown by thermal epitaxy (Fig. lb). (2) We have demonstrated that the electronic 
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Figure 1. Surface morphology of GaAs epitaxial layer 
electro-epitaxy; (b) by thermal epitaxy. 
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properties of electroepitaxial layers are superior to those grown, by standard 
methods. (3) The distribution coefficients of dopants (k Q ) a nd thus the 
carrier concentration can be readily controlled by the current density as 
shown in Fig. 2. (4) The growth velocity can also be controlled by the cur- 

rent density over a broad temperature range as shown in Fig. 3. In, addition, 
electroepitaxy lends itself to quantitative theoretical treatment; thus, in 
Figs. 2 and 3 the lines were theoretically obtained and are in good agreement 
with the experimental points indicated on the figure. (5) The dependence of 
the growth rate -on current density should permit the achievement of high growth 
rates and thus the growth of bulk GaAs material. (6) Since the current density 
controls the growth process, electroepitaxy, lends itself -uniquely to fundamental 
quantitative studies of crystal growth and segregation. 

The conclusion that bulk GaAs crystals might be obtained by electroepitaxy 
is currently under investigation. Our preliminary experiments showed that by 
electroepitaxy GaAs growth rates of about 20- ym/min can be achieved which are 
1-2 orders of magnitude greater than those achieved by standard LPE utilizing 
step cooling or the solution-substrate system. These growth rates are only 
about 1 order of magnitude smaller than those typically obtained in GaAs melt 
growth. It should be emphasized, however, that we have achieved the above 
growth rates employing our standard experimental system designed to .operate in. 
a rather limited range of temperatures and electric currents. Accordingly, in 
order to establish the real growth velocity limitations we have initiated the 
design, of an electroepitaxy system for growth under extreme conditions of electric 
current and temperature. Furthermore, we have initiated theoretical studies 
on the interface mechanisms limiting the growth velocity. This work is now in 
progress. We believe that electroeptaxy is a simpler and better suited method 
\for space processing of bulk GaAs crystals than melt growth, provided high 
growth rates -can be achieved. 
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Helt- Growth 

Growth from the melt is best suited for the growth of large single crystal 
within a relatively short period of time. However, as emphasized in our first 
annual report (April 1978) , extensive studies are required before device qualit 
GaAs can be obtained by this method. In this respect it should be emphasized 
that our recently developed characterization techniques have directly supported 
the conclusion that the quality of the commonly available GaAs is generally 

poor as manifested by the measured compositional and electronic inhomogeneities 

* * 

and by the high degree of compensation (see below). 

Presently our effort- is being focused on the initiation of melt growth 
experiments which will enable the investigation of the effects of As vapor 
pressure, the effects of ..convection in the melt, and the effects of growth 
interface phenomena on .characteristics .of GaAs single crystals. A melt-growth 
apparatus has been designed and is shown schematically in Fig. 4. It consists 
of a three-zone furnace designed to operate in a horizontal or in a vertical 
position. Construction of the apparatus has been completed in a "clean" labora- 
tory, and thermal characterization tests are now in progress. Following is a 
summary of equipment design, contruction and subsystem installation. 

Main Furnace Assembly : A complete set of plans were drawn and submitted 

to the Marshall Fiirnace Company for construction of the main furnace. There 
are several unique features incorporated in this furnace design: four indepen- 
dently controlled platinum heating coils, an internal cooling system utilizing 
an alumina tube to guide the cooling gas along the heated muffle; cooling gas. 
is supplied through a perforated inconel ring inside the furnace. This system 
generates increased thermal gradients without the risk of water-cooling tech- 
niques. High purity alumina tubing is used throughout the furnace to minimize 
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Support Systems 


Figure 4. Schematic representation ot the apparatus designed for the 
melt growth of GaAs ' 
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contamination and to provide fine adjustment of the thermal profile; power 
taps are placed in the three large windings. The furnace may be operated 
with or without standard size heat pipes. Installation <and initial checkout 
above 1000 °C have been completed. Stability and thermal profiling tests are 
forthcoming. 

Carriage, Track and Bench Assembly ; These systems provide physical 
support for the furnace, guidance during furnace movement and support for 
all of the furnace sub-assemblies. Design and contruction of all parts and 
sub-assemblies was carried out at M. I.T. Process tube mounts, drive cable 
guides, plumbing for gas and water, metering for process tube and cooling 
gases, thermocouple connection panels, furnace power connections and the 
cable drive unit were designed, built and installed. Plexiglas doors and 
shields are fitted .into rails on. the, bench .top, .then fastened .to the. fume hood 
to form a protective enclosure for the entire furnace assembly, low friction 
Thomson Roundway bearings support the furnace carriage on stainless steel ways. 
This method of guidance permits system operation in any orientation between 
horizontal and vertical position. Main support pivot, cable guide system and 
bench design permit rapid changes in system orientation with, only minor altera- 
tions. Key features of the support system are: isolation of the process tube 

through two layers of vibration damping composite material and an asbestos pad, 
adjustable process tube supports to compensate for various tube diameters, tube 
sag or other alignment variations and flexible high pressure tubing which allows 
cooling gas flow to the furnace under all operating configurations without modi- 
fications. A section of the bench top- may be removed to allow vertical operation. 
Adjustable indicators on the carriage are used to precisely position the furnace 


as well as .measure total furnace travel. 
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Gable Drive Unit ; This system was designed and constructed to operate 
under all anticipated loads generated by the furnace and carriage system. It 
consists of a D.C. control motor, controller, reductionigears and a threaded 
take-up drum. Stainless steel cable is used between the drum and carriage 
to minimize vibration transfer. The design of the cable guides permits the 
furnace track to be raised and lowered without changing the drive system or 
the furnace position. Speed in the drive motor is continuously variable from 
0.4 mm/.hr to 80 mm/hr in two overlapping ranges. Other speed ranges may be 
obtained by changing one gear pair. 

Feedthroughs, Thermocouples and Power Wiring ; Feedthrough endcaps were 
designed and built for two sizes of process tubes. These units are 0-ring 
sealed and water cooled (to protect the 0-ring). Overpressure of approximately 
0.5 atmosphere or moderately high vacuum may be sustained in the process tube 
using these endcaps. Each unit can pass up to seven thermocouple tubes, gas 
inlets/outlets or power cables. Quick release connectors and flexible tubing 
permits operational simplicity and vibration isolation from the input plumbing 
The rm o couple outputs are fed to two junction boxes, one at each- end of the 
furnace. Total capacity is 20 measuring thermocouples plus 4 control thermo-, 
couples. Both type K. (chromel-alumel) and type S (platinum-platinum-rhodium) 
thermocouples are used .for .the appropriate temperature ranges.- Power wiring is 
designed for simplicity of operation: no changes are required when altering 

furnace configuration and the use of polarized quick connect plugs eliminates 
wiring errors when disassembly is required. 

Monitoring .and Control Systems : Hard copy recording of thermal data- is 

provided by a 15 channel, time multiplexed chart recorder. This unit is 
linked directly to type K and type S thermocouple inputs and the buffered 



analog outputs from the four digital thermocouple readout displays* Precision 
millivolt sources were designed and built to offset each of the* analog outputs. 
This feature permits display of 1 to 50 C° full scale on the chart recorder. 
Furnace temperature is regulated by high sensitivity Leeds & Northrup con- 
trollers, driving SCR power packs. Each zone is independently controlled. 

Current metering and electronic current limiting plus a fail-safe feature 
in the control package prevent damage to the platinum furnace windings. 

Support Systems and Materials : A single zone furnace with control system, 

to be used for out gassing of source materials and ampoule while under vacuum, 

prior to melt-growth in the main furnace has been constructed. Vacuum is 

provided by a. cart mounted Varian diffusion pump system with ultimate pressure 
-7 

of 10 torr. This unit may be used for ampoule preparation or evacuation of 
the main furnace tube. Clean water for materials preparation is generated by 
the de-ionizing system. It consists of a carbon filter, demineralizer/ deionizer 
cartridge, heater (75 C°) storage tank, an activated carbon trap for organic 
impurities and two ion-exchange resin beds. Water of approximately 18 megohm 
quality is obtained from this unit. A sterile, 22 ym filter is installed at 
the outlet tap to trap particulates. 

Boats have been manufactured from high purity GE quartz. The boats are 
capable of holding approximately 35 g of material, with provision for a seed 
crystal, if desired. Two designs have been generated for the initial experiments. 
The ampoule has been modeled, with modifications in progress. Final dimensions 
will be determined after completion of temperature profiling tests. The ampoule 
is a two-piece quartz tube designed to allow independent evacuation of the source 
(arsenic) and boat (GaAs) regions at different temperatures. Fine diameter 
quartz wells protrude through the ampoule wall into the boat region. Thermocouples 
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placed in these wells will be used to monitor' thermal conditions in* the melt. 
(Direct contact of the melt and thermocouple is not permissible due to .the 
corrosive nature of arsenic at elevated temperatures.) Gallium, arsenic, 

-and gallium arsenide (polycrystalline) of greater tjian. 6N purity have been 
obtained from Metron, Inc. , with a spectrographic analysis of the GaAs boules 
included by the manufacturer. 

Growth-Property Relationships 

Very' rapid developments in our characterization techniques made pos- 
sible the initiation of direct experimental studies of growth-defect relation- 
ships earlier than originally scheduled. These studies were aimed towards ' 
the assessment of the role of the substrate as an active' source of point 
defects in epitaxially grown layers. 

Utilizing our SEM technique ‘we have measured the distribution of minority 

carrier lifetime in epitaxial layers, which provides a sensitive measure of the 

distribution of point defects acting as recombination centers. Typical results 

are given in Fig. 5; they clearly illustrate the importance of defect out- 

diffusion, from- the substrate into the ^epitaxial layer. Thus, the value of the 

lifetime in the vicinity of the epilayer-substrate interface is much lower 

(i.e. , defect concentration is much higher) than away from the interface. 

Furthermore, the lifetime decreases as the time during which the layers are kept 

at am elevated temperature increases. This, time dependence is consistent with 

the theoretical dependence describing out-diffusion of defects from the substrate, 

— 9 - *2 

and it was used to determine the value of defect diffusion constant, D - 5x10 an* / 

This work, initiated only recently, is still in progress. * The relationship 
between current density and carrier concentration (segregation) pointed out in 
elect roepitaxial growth above is also being pursued further. 



LIFETIME (s ) 


E( 

DIFFUSION 

CONSTANT 



o 5 . 


Figure 5. Lifetime as a function of time at growth temperature; crosses 
at substrate - epitaxial layer interface; circles away from 
interface. Lines correspond to out-diffusion of defects from 
substrate 
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CHARACTERIZATION OF GaAs 

As shorn in Table I, during the current year we have significantly 
advanced theoretical and experimental development of novel approaches^ to 
compositional , structural and electronic characterization of GaAs. We have 
also incorporated in our study standard advanced characterization methods. 

In the present report w e will highlight the unique aspects of our novel 
methods. 

Scanning IR Absorption 

Inhoraogeneities in semiconducting materials are commonly recognized as con- 
stituting a limiting major factor in microelectronic device applications. 
Nevertheless, the quantitative determination of inhomogeneities on a micro- 
scale have been so far possible only in elemental semiconductors (Ge and Si) • 

The most commonly used method to obtain' quantitative carrier distribution pro- 
files is based on spreading resistance measurements. Pressures metal contacts, 
required for this technique, are highly irreproducible in the case of compound 
semiconductors, and the technique is not applicable to GaAs characterization. 

Accordingly, we have developed a method and constructed an apparatus for 
the contactless and nondestructive determination of the spatial* variation of 
the free carrier concentration in GaAs - which is, of course, applicable to 
other semiconductors. The method is based on the quantitative relationship. 

between infrared absorption and free carrier concentration, and it is applicable 

15 —3 

to any semiconductor material with free carrier concentration exceeding 10 cm . 

The experimental arrangement of our system is shown in Fig. 6. 

A parallel beam from a >00^ laser (tunable between 9.16 and 11.02 pm)' is 

transmitted through the sample positioned on a stage with x-y motion and a 

13 

scanning rate ranging from 50 to 10 ym/min. The light transmitted through 



Table I 


rjetnoas Employed for GaAs Characterization 


Method 

Characterization' ■ 

Present Status 

Future Plans 

Novel Techniques 

- 

1 „ 


Scanning IR- 
Absorption 

Profiling of carrier concentra- 
tion & ionized impurities 

Apparatus completed, theory devel- 
oped, results obtained for n-type 
GaAs 

p-type GaAs 

SEM-Electron Beam 
Induced Current 

Instantaneous profiling of dif- 
fusion length, lifetime & sur- 
face recombination velocity 

Apparatus completed, method 
developed & demonstrated 

Quantitative analysis, 
experimental extension to 
low & high temperatures 

Derivative Surface 
Photovoltage Spec- 
troscopy 

Determination of optically ac-’- 
tive defect impurities &.band 
structure transitions 

Method developed, preliminary 
Results obtained 

Transient analysis yielding 
dynamic parameters of deep 
centers; extension to inter- 
face characterization 

Transport Tech- 
nique 

Determination of compensation 
ratio & total amount of 
ionized impurities 

Theory developed for n-type GaAs, 
results obtained 

Theoretical extension to 
minority carriers - in p-type 
GaAs 

Standard Techniques 




Cathodoluminescence 

Spectroscopy 

Profiling of defects & impuri- 
ties active In radiative re- 
combination 

Apparatus components obtained, 
assembly in progress 

Cross-correlation with 
EBIC technique 

Photoluminescence 

Spectroscopy 

Precise determination of radia- 
tive transitions, identifica- 
tion of residual impurities & 
defect complexes 

Experimental arrangement designed, 
apparatus components being 
acquired 

Cross-correlation with deri-' 
vative surface photovoltage 
spectroscopy & TDLS 

Transient Deep 
Level Spectros- 
copy (TDLS) 

Determination of energy levels 
& dynamic parameters of deep - 
centers 

Apparatus components acquired, 
assembly in progress 

Cross-correlation with deri- 
vative surface photovoltage 
spectroscopy 




Figure 6. Schematic representation of the apparatus.- for IR-scanning absorption 
spectroscopy 
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the desired area of the sample is focused on. a detector (bolometer) through 
the optics of an IR microscope; the signal from the detector, which is pro- 
portional to. the intensity of radiation, is amplified and recorded. 

The system permits a continuous monitoring of transmittance of the small 
wafer area, which is converted, through a data processing system, into 
variations of carrier concentration and of the compensation ratio. We have 
recently completed a rigorous theoretical treatment of the pertinent optical 
phenomena in n-type GaAs required for the quantitative processing of the data. 
A typical microprofiling of GaAs obtained recently is given in Fig.. 7. It is 
seen that a two-dimensional quantitative profile of the carrier concentration 
in GaAs single crystals.- has been achieved for the first time. This profile 
clearly shows the pronounced inhomogeneities in the presently available melt- 
grown GaAs. 

The spatial resoluton of our IR scanning system is about 20 ym and varia- 

*t / n 

tions in the free carrier density in GaAs as low as 5 x 10 cm" in n-type 

14-3 - — 

and 2.5 x 10 cm in .p- type material can be detected. 

Employing the same experimental system and a theoretical model wer.have 
developed we can determine in GaAs the average value of the compensation ratio 
as well as its variations. 

SEM Techniques 

In our approach scanning electron microscopy (SEM) is used primarily 
in the electron beam-induced current (EBIC) mode. Thus, the electron beam 
represents in essence a highly focused excitation source for excess minority 
carriers which are collected by a Schottky barrier or a p— n junction. 

We have installed an SEM and accessory electronics especially designed 
for our EBIC mode approach (see Fig. 8). In parallel, we have been pursuing 
a theoretical and experimental study to advance our understanding of the 
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Figure 7. Typical carrier concentration microprof ile of melt-grown 'GaAs; 
results obtained with IR-scanning absorption technique 
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excitation and recombination processes associated with the EBIC mode. In 
these studies we have used Si as a convenient model material, and we have 
obtained and confirmed high resolution profiles such as the minority carrier 
diffusion length, the lifetime, and in certain instances, the- dopant profile. 
Tupical results illustrating these possibilities are shown in Fig. 8 together 
with a schematic arrangement of the experimental system. For the first time 
a direct correlation has been achieved on a microscale between compositional 
inhomogeneities- and electronics parameters such as minority carrier- diffusion 
length and lifetime. 

We have .utilized this technique to analyze the lifetime distribution in 
epitaxially-grown GaAs. The results (shown earlier in Fig. 5) provided direct 
evidence of out-diffusion of defects from substrate to epitaxially grown layers. 

Extensive experimental work is in progress devoted to the extension of 
EBIC-mode technique to temperatures well below and well above room temperatures. 
In spite of inherent technical difficulties, we consider such an extension of 
key importance in identifying the specific mechanisms of minority carrier recom- 
bination which are the controlling factors in device performance. We have also 
acquired apparatus and accessories for carrying out cathodoluminescence spec- 
troscopy employing our SEM system. The assembly of this apparatus is in 
progress . 

Derivative Surface Photovoltage Spectroscopy 

Most recently we have developed a derivative form of Surface Photovoltage 
Spectroscopy based on the wavelength, modulation of the incident light. In its 
standard (nonderivative) form surface photovoltage spectroscopy was successfully 
utilized, however, primarily for surface state studies. 

Derivative photovoltage spectroscopy makes possible, in a single experi- 
ment, the quantitative characterization of optically active deep levels 
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(impurities and defects) , of interface states, and of the energy band structure 
above the energy gap. Preliminary results are shown in Figures 9 and 10. Figure 
9. shows deep level transitions associated with impurity traps, whereas Fig. 10 
shows critical point transitions. It should be emphasized that there exists 
hardly any other technique permitting the simultaneous characterization of 
these transitions. 

In conjunction with these findings we have developed a data processing 
system (based on a theoretical quantum defect model of photoionization transi- 
tions) which permits a direct determination of such defect parameters as energy 
position and relative concentrations from the measured first-derivative subband- 
gap spectra. 

Extensive studies on this technique are in progress and more detailed 
results will be presented in our annual report.' The same experimental system', 
employed in derivative surface spectroscopy will be complemented with apparatus 
for simultaneous photoluminescence and Transient Deep Level Spectroscopy measure- 
ments.. The experimental arrangement has already been designed and the apparatus 
components are being acquired. 

Transport Technique 

We have completed the development of a characterization approach based 
on low- and -high-frequency transport phenomena in n-type GaAs. Thus, utilizing 
an- advanced theoretical treatment in conjunction with experimental studies of 
the electron mobility and of the free carrier absorption we succeeded in 
developing reliable, practical procedures for the determination of the compen- 
sation ratio and of the total concentration of ionized impurities in GaAs. 

We utilized this approach in order to evaluate the quality of available 
melt-grown GaAs. The results are shown in Fig. 11. The 0=0 line (0 is the ■ 








2 (orb. units) 
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Figure 10. Band structure critical point transitions in GaAs, as revealed by 
second derivative surface photovoltage spectroscopy 




. Figure 11. Election mobility of GaAs- as a function of electron concentration 
calculated for different values of compensation ratio. Shaded 
area refers to the mobility range of commercial GaAs. 
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compensation ratio) represents the theoretical mobility limit. The shaded 
area represents the mobility range of presently available bulk GaAs. It is 
evidence that GaAs performs below the theoretical limit and that this tendency 
is particularly pronounced for higher purity (lower electron concentration) 
material. On the basis of our theoretical and experimental results we have 
developed most useful tables from which the compensation ratio can be directly 
obtained from carrier concentration and from mobility or absorption coefficient 
measurements. These tables are far more accurate than those in use up to this 
time. 

It has been recently reported that the compensating centers in GaAs are 
associated with point defects and defect- impurity complexes. Thus, our method 
for the determination of the compensation ratio will be used as a convenient 
means for relative determination of the overall defect density in GaAs. -We 
plan to utilize this method in all stages of our program in direct correlation 
with growth experiments. 



Electron mobility and- tree-carrier absorption in GaAs: 
Detetminatibh of the compensation ratio 
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Theoretical calculations of electron mobility and free-carrier absorption in n-type GaAs 
at room temperature were carried out taking into consideration all major scattering 
processes. It was found that satisfactory agreement between theoretical and 
experimental results on free-carrier absorption is obtained only -when the effect of 
compensation is quantitatively taken into account. In conjunction with experimental 
studies it is shown that the electron mobility (for n > 10 1 5 cm -3 ) and , free-carrier 
absorption (for n > 10 16 cm” 3 ) are sufficiently sensitive to the ionized impurity 
concentration to provide a reliable means for determining the compensation ratio. 
Convenient procedures are presented for the determination of the compensation ratio 
from the free-carrier absorption coefficient and from the computed values of room- 
temperature electron mobility. Values of the compensation ratio obtained by these two 
procedures are in good agreement provided the carrier-concentration variations in the 
material are not appreciably greater than 10%. 

\ 

PACS numbers: 72.20.Dp, 78.20.Dj, 72.20.Fr, 72.80.Ey 
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I. INTRODUCTION 

III-V semiconductors usually .exhibit an appreciable 
degree of compensation. Knowledge of the extent-of com- 
pensation is important for the understanding of the basic 
- electronic characteristics of these materials and their devices 
and.for the .evaluation of crystal-growth processes. The 
methods for determining the compensation are usually relat- 
ed to the fact that the compensation reduces the free-carrier 
concentration and it enhances free-carrier scattering by ion- 
ized impurities. 

In GaAs, because of the small donor ionization energy 
5'meV), the accurate determination of donor and accep- 
tor concentrations from free-carrier concentration changes 
must be carried out at temperatures well below 15 °K. Simi- 
larly, the determination of the compensation ratio from the 
amplitude of Schubnikov-de Haas oscillations, as proposed 
recently, 1 must be carried out at 4.2 °K. The low-tempera- 
ture requirements render these methods inaccessible for rou- 
tine needs. Instead, methods basedon scattering of free ear- 
ners by ionized impurities are employed. The reliable 
determination of the compensation ratio from' free-carrier 
scatteringby ionized impurities is complicated', however, by- 
the fact that longitudinal optical phonons contribute signifi- 
cantly to the overall scattering. 2 For this polar-mode scatter- 
ing, the relaxation time cannot be defined except for low 
temperatures.* Accordingly, it is necessary.to use variation- 
al 4 or iterative procedures 2 to combine all relevant scattering 
mechanisms. Earlier approaches 5 6 to the determination of 


°On leave from Institute of Physics. Polish Academy of Sciences, Warsaw. 
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OF POOR QUALITY 

the compensation ratio from electron mobility neglected 
these aspects. More recently, a rigorous semiempirical pro- 
cedure has been formulated 7,8 for determining the compen- 
sation ratio in /z-type GaAs from the Hall constant and resis- 
tivity measurements at 77 °C. This procedure.takes into 
consideration all major scattering mechanisms-but it is ap- 
plicable only to nondegenerate material. Theoretical calcu- 
lations of electron mobility (based on an iterative procedure 
for solving the Boltzmann equation) including compensa- 
tion as a parameter and considering all major scattering 
mechanisms have also been reported. 2,9 Although the screen- 
ing of polar vibrations by free carriers (which is significant at 
room temperature for n > 10 !7 cm' 3 ) was not taken into con- 
sideration, these results represent the most reliable theoreti- 
cal values of electron mobility in GaAs thus far. 

Free-carrier absorption has been used also for the deter- 
mination of the compensation ratio in /z-type GaAs. 10 " 12 In 
this approach (in contrast to that based on the electron mo- 
bility) the ionized impuritycontnbution-to free-carrier ab- 
sorption is v separable from contributions by other scattering 
mechanisms. This approach has not as yet been correlated 
with compensation-ratio values obtained by other methods. 
Furthermore, the values of free-carrier absorption available 
in the literature are not very reliable because the appropriate 
material parameters were not used in carrying out the 
computations. 1314 

In the present study; the effect of compensation on two 
independently measured quantities, i.e., on the electron mo- 
bility and on the free-carrier absorption, is considered Theo- 
reticahcomputations of mobility are based on a variational 
method in the form proposed in Ref. 4 and they include all 
major scattering processes and.. screening effects Computa- 
tions of the free-carrier absorption coefficient arc based on 
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the approach of Ref. 10. A comparison is made with experi- 
mental results and with data available in the literature 
A procedure is outlined for the convenient determination of 
the compensation ratio from measured values of room-tem- 
perature electron mobility and from the absorption coeffi- 
cient at a wavelength of about 10 /zm. x 

II. THEORY 

A. Genera! considerations 

The present theoretical calculations are intended to^ 
provide the basis for the determination of the compensation 
ratio (or total density of ionized impurities) in n ~ type GaAs 
from independent room-temperature measurements of the 
free-carrier absorption and the electron mobility. 

The electron mobility /z and.the free-carrier absorption 
coefficient a will be defined as 

{i=cr(0)/Tie, a{co) — A7Tu{cd)/cn ^ (1) 

where cr(0) and a{cci) are’the dc and ac conductivity, respec- 
tively, n is the concentration of free electrons, c is the.velocity 
of light, and n r is the refractive index at a given irradiation 
frequency oo. 

Whemthe energy-band structure, the electron-phonon 
coupling constants, and' the electron concentration are 
known, the value of cr can be calculated if the total concen- 
tration of ionized impurities is also known. For //-type 
GaAs, the available material parameters are sufficiently ac- 
curate for a quantitative calculation of / 1 and a{co). For a 
given temperature, such calculations provide two absolute 
quantities, i.e., /z and consf , and one functional depen- 
dence a{o)\ each of whichcan be used to determine the con- 
centration of ionized impurities. Accordingly, in cases 
where ionized impurity scattering is appreciable, experimen- 
tally, measured electron mobilities, free-carrier absorption 
coefficients, and the frequency dependence of absorption co- 
efficient can serve as the basis for testing the validity of the 
theoretical calculations and, the reliability of the procedure 
for'the determination of the compensation ratio. 

In general, dc and ac conductivities can be treated in the 
same theoretical framework. In practice, however, it has 
been shown that dc conductivity is satisfactorily obtained 
from the solution of the semiclassical Boltzmann equa- 
tion, 2,4,19 whereas ac conductivity can be conveniently ob- 
tained from second-order perturbation theory. 10 - 20 - 2 ? 'Both 
approaches are utilized imthe present study; the computa- 
tions are carried out using the material parameters given in 
Table I and discussed in detail in Ref. 2. 

In view of the high (1.43 eV) and directienergy gap of 
GaAs, a single spherical and parabolic conduction band de- 
scribed by an effective mass m* will be assumed. With such 
an approximation, the electron mobility is slightly overesti- 
mated and the free-carrier absorption is underestimated. For 
this reason an effective mass m*/m 0 = 0.068 is used, which 
is somewhat higher than the exact value of the effective mass 
= 0.064 ± 0.002 at the band edge of GaAs. 23 

’ * In the theoretical approach adopted in' the present 
study, dilute solid solutions are assumed; i.e,, interactions 
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TABLE I GaAs parameters used in present computations (300 *K) {after 
Refs 2 and 9) 


Low-Uequency 
dielectric constant c n 
High-frequency 
dielectric constant e. 
Optical-phonon 
energy fuj 0 
Deforma lion 
potential E t 
Longitudinal 
elastic constant C, = pv) 
Piezoelectric 

coefficient h ~ U {5/Cj+ 4 /C f ) 
Effective mass m*/m 0 


between impurities, quantum corrections for ionized impu- 
. rity scattering, 24 and' differences in the short-range part of 
the impuritypotential are not taken into consideration. This 
assumption sets a limit to the electron (and impurity) con- 
centration, which for / 2 -type GaAs at room temperature can 
be realistically estimated to be 3 X 10 l * cm' 3 . In this respect, it 
should be pointed out that 1 differences in scattering by var- 
ious impurities (Te, Se, and S) in GaAs at concentrations 
above 10 IS cm" 3 have been indicated in Refs. 25 and 16-18. 
However, a study of the electron mobility in about 10 3 GaAs 
single crystals doped with Te, Se, Sn, and Si has failed to 
show any consistent differences among the various do- 
pants 26 ; the only difference (slightly lower mobility) was 
found in Si-doped crystal (at concentrations exceeding 10 18 
cm" 3 ); this result was apparently due to the amphoteric be- 
havior of Si. 

B. Free-carrier absorption 

The present analysis will be confined to wavelengths 
greater than 4 /zm, where transitions between various con- 
duction-band minima are negligible. Thus, only the lowest 
conduction band with a minimum at the center of the Bril- 
louin zone (T 6 ) needs to be considered. The absorption of. 
free carriers in a single conduction band constitutes an indi- 
rect transition in which momentum is conserved through 
interactions with the lattice. Using the approximation of a 
parabolic conduction band, the absorption coefficient can be 
directly calculated from the expressions derived in Ref. 10. 

The total absorption coefficient a f is. 

1a,=a op +a at +a, mp , - - (2) 

where a opi . a ac , and a imp are the absorption coefficients 
corresponding to electron interactions with screened optical 
phonons, acoustic phonons, and screened ionized impuri- 
ties, respectively. Piezoelectric scattering is not taken into 
consideration; mobility calculations (see below) showed that 
for /z-type GaAs, at room temperature thecontribution of 
piezoelectric scattering to the total mobility is less than 2%; 
a similar contribution (which is within the accuracy of the 
calculations) is expected to the free-carrier absorption. 

According to Ref. 10, the absorption coefficients a v can 
be written as follows,, where v designates op, ac, or imp: 

a v - E f [A*)-A x' v )]dx, (3) 

z 3 Jo 
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1291 
10 91 
36 meV 
7 eV 

14.03X10" dvn/cnr 

2 39 X 10-- 
0.068 (see tevt) 
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where J{x) = l/[exp(x— 77)+ 1] is the electron distribution 
function x—E/k 0 T, y=Ej/k<J\ and z=-iuo/k Q T arc the re- 
duced.electron energy, the Fermi energy, and the incident- 
photon energy, respectively, k* is the Boltzmann constant*, 
and T is the absolute* temperature The values of' 
C v , G yy and x v are for optical phonons 


C OR =2.96Xl0 s (mVm 0 ) 


n r € 0 J 


X 


sinh-z 

2 


sinh^z, sinh-z' 
G ox ~2P op [\+ ~)-aL 


x op =x+z±z„ 
for ipriized impurities, 


N, 


imp 


C imp =2. 175X10-’— — , 


^irrip 4 ^imp & 


x imp =x+z. 


imp. 


(4) 


(5) 


and for acoustic phonos. 


CL = 4.79 X 10 1 


,/jw* y te\ 

\ iri 0 ) n r pvf 


G ac ~ QRJGlWAL PA.GE SS 

x ac =x + z, oF POOR QUALITY 

where 

P v — (xx \) 1/2 , B v = x+ x\„ 

B v + 0 + 2 P v 


.= ln(- 


Bf - 2P V 
■D v =(B v + ay-4Pl 


)■ 


( 6 ) 


m 


The symbols in Eqs. (4)-(.7) are as follows: € 0 and are the 
static and the high-frequency dielectric constants, respec- 
tively, N imp = N £ -}- iVj is the concentration of ionized 
impurities, E x is the conduction-band deformation potential, 
p is the crystal density, v l is the longitudinal sound velocity, 
z { — ficot/koT is the reduced energy of longitudinal optical 
phonons, m 0 is the electron mass, a is the reduced screening 
energy, 

a = ft/2m*iyc 0 T, 

where the screening length I D is defined as 

= 5.80x 10” 

ll 


( 3 ) 


^0 


(9) 


and F t \(y) is the n - th order Fermi Dirac integral 
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It* shouldbe noted that a op is a sum of t\vo«terms given 
by Eq ’(3) which correspond to x op = x 4- z — z t and ^ 
z* — z — z 7 andjc^ p = x 4 - z 4- z /( andz' — z 4- z b respectively. 
A discussion of the numerical calculations carried out on the 
basis of Eqs. (2)-(9) and the GaAs parameters listed in Table 
I will be presented below together with experimental results. 

The carrier concentration n — N £ — N J* (where 
N p and N are the concentrations of ionized donors and 
acceptors, respectively) enters into Eqs.-(2)-(9) througluthe 
reduced Fermi level 77 according to 


«=(2/7T ,/2 )yV c F l/2 ( 7 /) 


(10a) 


or 

n = 5.44 X 10 15 [ 7’(ot */ m 0 ) ] 3/2 F y/1 (jf), - (10b) 

where N c is the effective density of states in the conduction 
band. 


The compensation ratio 6 is defined as 6—N^YN 
accordingly, the total concentration of the ionized impuri- 
ties is related to the compensation ratio as follows: 


N^-nd+eyo-B). . (H) 

It is apparent from Eqs. (2)-(9) that for a given frequen- 
cy of the incident light and for a given electron concentra- 
tion, only the term a imp depends on the compensation ratio 
(being proportional to A^ mp ). thus, if a eNp is the measured 
absorption coefficient in a material with a compensation ra- 
tio 6y and <2 0p , a ac , and a; mp are theabsorption coefficients, 
computed [from Eqs. (3)-(9)] for the same wavelength and 
the carrier concentration taking 6—0 (i.e., iVj mp =w), the to- 
tal concentration of ionized impurities can be determined by 




exp (^op“F^ac) 




imp 


(12a) 


which according to Eq. (1 1) leads to the following expression 
for the compensation-.ratio: 


q a t\ p ( g opri’^ac^F g imp) (pb) 

^*cxp "F^imp fa’op^^ac) 

The values of a imp , a opy and a ac computed as a function of 
n % for a wavelength of 10 are, given in Table II. 

A direct comparison between the experimental values 
of absorption coefficient and the theoretical dependence of 
the total free-carrier absorption [computed from Eqs. (2)- 
(7)] on electron concentration is not possible since samples 
withfvarying electron concentration can be characterized by 
a different degree of compensation. To provider means for 
such.comparison, a reduced experimental absorption coeffi- 
cient, a* xp will be introduced representing a exp for zero 
compensation. 


• — 1 ®cxp (13a) 

which according to Eq. (l l) becomes 

The theoretically obtained data in Table II can be used for 
relating a* xp to the compensation ratio as determined by oth- 
er means (e.g., by electron-mobility measurements). This 
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TABLE II Compuied values of electron' mobdily and ir absorption in /My pc GaAs 


Compen- 
sation Electron mobiHtv (cm : /V sec) Absorption 

ratio * at 10 urn (cm' 1 ) 


0 = 

00 

0.1 

02 

03 

0.4 

05 

0.6 

0.7 

08 

09 

®inip 


cr 

1 X io” 

7810 

7760 

7700 

7620 

7530 

7400 

7220 

6940 

6490 

5590 




1.5X10” 

7740 

7680 

7590 

7490 

7360 

7200 

6970 

6640 

6120 

5150 




2X10” 

7600 

7510 

7410 

7290 

7140 

6950 

6690 

6330 

5770 

47S0 




3X10” 

7540' 

7430 

7300 

7150 

6960 

6730 

6420 

6000 

5400 

4370 




4X10” 

7450 

7320 

7170 

6990 

6780 

6510 

6 ISO 

5730 

5100 

4070 




5X10” 

7370 

7230 

7060' 

6S70 

6640 

6350 

6000 

5540 

4850 

3870 




6X10” 

7280 

7130 

6950 . 

6730 

6490 

6190 

5820 

5340 

4700 

3670 




7X10” 

7230 

7060 

6860 

6640 

6380 

6060 

5680 

5200 

4540 

3510 




8X10” 

7150 

6970 

6760 

6530 

6260 

5930 

5540 

5050 

4390 

3370 




9X10” 

7050 

6860 

6650 

6400 

6120 

5790 

5390 

4900 

4250 

3232 




1 X 10” 

6980 

6790 

6560 

6310 

6020 

5680 

5280 

4790 

4140 

3120 

0 0036 

0 025 

0.341 

1.5X10” 

6710 

6*80 

6230 

5950 

5640 

5280 

4870 

4370 

3720 

2710 

oooso 

0 038 

0.510 

2X10” 

6500 

6250 

5980 

5680 

5360 

5000 

4580 

40S0 

3430 

2440 

0.0142 

0 050 

0.679 

~ 3X10” 

6190 

5920 

5630 

5320 

4980 

4610 

4180 

3680 

.3030 

2060 

0.032 

0 075 

1.10 

T £ 4x10” 

5970 

5690 

5380 

5060 

4720 

4340 

3910 , 

3410 

* 2760 

1820 

0.056 

0 100 

1.35 

~ 5X10” 

5810 

5510 

5200 

4S70 

4520 

4140 

3710 

3200 * 

2560 

1640 

0.087 

0 126 

1.68 

i 6X10” 

5680 

' 5370 

5050 

4720 

4370 

3980 

3550 

3040 

2400 

1510 

0.125 

0 152 

2.00 

8 7X10” 

5570 

5260 

4940 

4600 

4250 

3860 

3420 

2910 

2270 

1400 

0.169 

0.178 

2.32 

g 8X10” 

5480 

5160 

4840 

4490 

4140 

3740 

3300 

2790 

2160 

1320 

0.220 

0.203 

2 65 

| 9X10” 

5400 

5080 

4750 

4400 

4040 

3640 

3200 

2690 

2060 

1240 

0 277 

0.229 

2.97 

w lxio” 

5330 

5010 

4670 

4320 

3960 

3560 

3120 

2600 

1980 

1180 

0 340 

0.255 

3 28 

1.5X10” 

5090 

4750 

4400 

4030 

3650 

3240 

2790. 

2280 

1690 

960 

0 746 

J&387 

4 84 

2X10” 

4910 

4540 

4160 

3780 . 

3390 

2970 

2520 

2030 

1470 

820 

1.29 

Ck52I 

6 35 

3X10” 

4730 

4350 

3960 

3560 

3160 

2730 

2280 

1790 

1270 

680 

2.79 

0.796 

9.26 

4X 10” 

4530 

4130 

3730 

3320 

2920 

2490 

2050 

1590 

1100 1 

580 

4.82 

1.08 

12.1 

5X 10” 

4470- 

4060 

3640 

3220 

2800 

2380 

1940 

• 1490 

1020 

530 

7.26 

1.34 

14.8 

6X10” 

4370 

3950 

3520 

3100 

2680 

2260 

1830 

1390 

950 

480 

10.2 

1.68 

17.5 

7x10” 

4290 

3860 

3420 

3000 

2580 

2160 * 

1740 

1320 

■ 890 

450 

„ 13.5 . 

1.99 . 

20.1 

8x10” 

4220 

3780 

3340 

2910 

2500 

2080 

1670 

1260 

. 840 

420 

17.2 

2.32 

22.7 

*9X10” 1 

4160 

3710 

3270 

2840 

2430 

2010 

1610 

1200 

800 

400'" 

21.3 

2.65* ' 

25.3 

IX 10” 

4100 

3640 

3200 

2770 

2360 

1950' 

1550 

1160 

770 

380 

' 25.8 

3.00 

27.9 

1.5X10” 

3860 

3400 

2950 - 

2520 

2130 

1740 

1370 

1010 

660 

330 

53.4 

4.85 

40.5 

* 2X10” 

3690 

3220 

2780 

2360 

1980 

1600 

1250 

920 

600 

290 

88.8 

6.91 

52.8 

3x10” 

3460 

3000 

2560 

2160 

1790 

1440 

1120 

810 

530 

260 

178 

11.6 

76 6 

4X10” 

3300 ‘ 

2840 

2410 

2020 

1670 

1340 

1040 

750 

480 

230 

284 

16.9 _ _ 

98 6 

5X10” 

3200 

2750 

2330 

1950 

1600 

1280 

990 

710 

460 

220 * 

400 

22.8 

120 


procedure will be discussed below in conjunction with ex- 
perimental results. At this time, an analytical extension of 
the above procedure to other wavelengths will be considered 
without resorting to additional numerical calculations. In 
previous treatments / 11, 20 ~ 2 \ } approximate analytical expres- 
sions have been obtained on the frequency dependence of the 
absorption coefficient assuming a nondegenerate electron 
gas. They were found to be in fairly good agreement with 
experiment even in the case of highly doped materials for 
which the assumption of nondegeneracy is not satisfied. It 
can be shown, however, that these results are not inconsis- 
tent since similar analytical expressions can be, obtained 
from Eqs. (2)-{9). 

Thus, in a spectral region, where fia is much greater 
than fia) Q , Ep, k<>T> and ak 0 T , Eqs. (3)-{7) reduce to 


a 0 ~-^fF l/2 C> 7 ) 

z 



These expressions give the same frequency (or wavelength) 
dependence of the absorption coefficient as the expressions 
in Ref. 1 1 and 20-22. For GaAs at room temperature andin 
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the spectral region of4— tO/zm (inequality fuo^fuo^ fc n k u i , 
and ak^T is fairly well satisfied up to an electron concentra- 
tion of about 10 1H cm - ’. For higher electron concentration 
Eqs. (14a)-(14b) still constitute a reasonable approximation 
for shorter wavelengths (e.g., 4-7 /zm), in agreement with 
previously reported experimental results. 111 ' 

According to Eqs (14a)-(l4c), the procedure for deter- 
mining the compensation ratio, from Eq. (12), and the nu- 
merical values m Table II for A o =10/zm, can be simply ex- 
tended to other wavelengths A by calculating the absorption 
coefficient according to the approximate relationships 

aoM. =a op^ l/A o) 25 ' (15) 

For an electron concentration n — 10 iS cm -3 and for 8<A< 12 
ftm the error in absorption coefficients calculated according 
to Eq. (15) [rather than from more exact numerical calcula- 
tion based on Eqs. (3)— (7)] does not exceed 9% for 
a imp and 4% for a op +a ac . These errors decrease rapidly 
with decreasing electron concentration. 


C. Electron mobility 


The present computations of electron mobility are 
based on a variational principle method in the form pro- 
posed in Ref. 4 and 19. This method allows the combination 
of all scattering mechanisms without invoking Matthiessen’s 
rule, and it has been successfully applied to the calculation of 
mobilities in a number of semiconductor compounds. 319 27-30 

The following scattering mechanisms are considered: 
screened optical-phonon scattering; 31 screened iohized-im- 
punty scattering, 19 * 32 acoustic-phonon scattering-through 
the deformation potential, 33 and piezoelectric interactions. 34 
The final expression used (see the Appendix) is 



(16) 



FIG. 1 Comparison of elivlP'm moKhties calculated by a va national meth- 
od (present calculation'^ and an iterative procedure ( Ref. 21 for /My pc 
GaAs at room icmp--rati*rc. ncreenmg of polar scattering ji» not included 
(seeieM) 



Electron concentration (cm -3 ) 


FIG. 2. Calculated component and total electron mobility as a function of 
electron concentration in «-type GaAs at room temperature. 


where D y/ 1%in and D are the determinants defined in the 
Appendix. 

Mobility values as a function of electron concentration 
obtained in the present study are shown in Fig. 1 together 
with results reported in Ref. 2 where an iterative procedure 
was used for solving the Boltzmann equation; screening of 
polar scattering intentionally is not included in the calcula- 
tions shown in Fig. 1 since this type of screening was not 
included in the calculations of Ref. 2. It is seen that the two 
methods are in very good agreement; the observed discrep- 
ancy at low and high electron concentration is primarily due 
to the assumption of a parabolic conduction band used in the 
present study; in Ref. 2 a nonparabolic conduction band was 
used with a slightly smaller band-edge effective mass 
(0.066). 

The component electron mobilities [i.e., electron mobil- 
ities calculated from Eq. (16) for each scattering mechanism 
acting alone] v are shown as a function of electron concentra- 
tion-inFig. 2. It is seen that the maximum contribution to the 
total mobility from piezoelectric scattering does not exceed 
2%, consistent with the assumption made above in the case 
of free-carrier absorption. It is also seen that screening ef- 
fects become pronounced for an electron concentration ex- 
ceeding 10 17 cm“ 3 . 

In Fig. 3 the total electron mobility calculated from Eq. 
(16) is shown together with the electron mobility obtained by 
summarizing the component mobilities according ro Matth- 
iessen’s rule(i.e M 1/// — 2,1/zq) It is seen that for an electron 
concentration of 10 17 cm - ' Mattheisscn's rule leads to an er- 
ror of about 30 %. Thus, in determining the compensation 
ratio one.must reiy on the results of numerical < ..cuhuieiis 
of mobility as a function of earner concentration iuc v n- 
pensarion ratio. The results obtained in the present ut , re 
summanzed in Table II. 
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FIG 3. Comparison of electron mobilities in /i-type GaAs calculated by a 
variational method and by Matthiessen's role employing the component 
mobilities. 



FIG 4 Theoretical (solid lines) and experimental values of mobiliiv as a 
function.of electron concentration m w-type GaAs for various compensa- 
tion ratios: the mobility range for commercially available GaAs is also 
indicated. 


III. COMPARISON WITH EXPERIMENT 
A. Experimental procedure and sample 
homogeneity 


mobility can be higher 36 and the infrared absorption lower 
than theoretical limits. Accordingly, only samples with car- 
rier-concentration inhomogeneities not exceeding approxi- 
mately 10% are considered in this study. The parameters of 
these samples are listed in Table III. 


Hall effect, dc conductivity, and infrared absorption 
measurements were carried out on /i-type GaAs samples 
with electron concentrations ranging from 6x 10 15 to 
2.6 X 10 18 cm' 3 . 

Melt-grown GaAs single crystals obtained from com- 
mercial suppliers and two epitaxially grown Te-doped layers 
(intentionally compensated with Ge) were used in this study. 
In order to obtain consistent and meaningful results, it was 
found necessary to assess the homogeneity of the samples. 
The carrier concentration distribution on a macroscale and 
microscale was determined employing a recently reported ir 
scanning technique. 35 In the melt-grown commercial crys- 
tals, the spatial variation of the absorption coefficient at 
/i= 10.6/mi was found to be of the order of its average value 


The experimental values of the electron mobility and 
electron concentration were obtained, from Hall-effect and 
conductivity measurements. In order to avoid ambiguities 
associated with the value of the Hail factor r the Hall con- 
stant was determined as a function of magnetic field up to a 
value of 100 kG; the saturation value (which corresponds to 
r— 1) was used to determine the electron concentration. 

The results obtained are shown in Fig. 4 together with 
the theoretical dependence of electron mobility on electron 
concentration for different compensation ratios. The results 
of Ref. 18 which will be discussed in conjunction with free- 
carrier absorption are also shown. It is seen that all mobility 
values are well below the theoretical values at zero compen- 
sation. From the electron concentration and the electron 


(or An of the order of n t 

3U ). In such samples the measured 

mobility, the compensation ratio for each sample was deter- 

r 

TABLE III. Compensation ratios derived from mobility and ir absorption 

measurements 

• - 


Sample 

Electron 



Homogeneity 

Mobility 

Compensation ratio 9 

No. 

Cone. (cm’ J ) 

Growth 

Dopant 

An/n ^ 

(cnr/V sec) 

Mobility 

Absorption 

1 

3X10“ 

Melt 

Te 

< 10 % 

2700 

0 S3 + 0.03 

. 0.89 - 0.05 

2 

1 . 1 X 10 ” 

Melt 

7 

< 10 % 

. 1800 

0.83 4- 0.03 

0.85 + 0 03 

3 

5.5X10” 

Melt 

Te 

< 10 % 

3300 

0.25 + 0 06 

0.23 + 0.08 

Z \ 

8.5X10" 

LPE 

Te 

<3% 

1350 

0.67 ± 0.03 

0 65 + 0 04 

5 

1.5X 10“ 

LPE 

Te 

<3% 

800 

0.76 + 0 02 

0.77 + 003 

6 

2 . 0 X 10 “ 

Melt 

7 

<15% 

2100 

0.37 + 0.08 

0.29 + 0.08 

7 

6 X 10 " 

Melt 


< 20 % 

3300 

>0.9 

>09 

8 

2.6 X !0 U 

Melt 

... 

< 20 % 

1650 

0.40 ±0.1 

0.3 0 I 

Ref. 18 

6 X 10 " 

... 

Te 

... 

3600 

0.18 

0.18 

Ref. 18 

1.15X10" 

... 

S 


2600 

0.31 

0.25 

Ref. 18 

1 . 2 X 10 " 


Tc 

... 

3100 

0.22 

0.22 


^Intentionally compensated with Ge (see text). 
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mined directly from Table II, the values of 0 are given in 
Tabic III. The normal range of mobility m commercially 
available high-quality GaAs single crystals us indicated in 
Fig 4 by the striped area; in the low-electron-concentration 
region, the available GaAs is signifkciUly compensated and 
typically the total concentration of iom::ed impurities ex- 
ceeds by an order of magnitude the electron concentration. 

The absorption coefficient was determined from the 
transmittance T (measured with a Fourier spectrometer) 
using 

T _ (l-i?) : exp(-Qt/) - 17 . 

1— ^? 2 exp(— 2acf) ’ V ' 

where d is the thickness of the sample arid R is the reflectance 
of GaAs (e.g., forA= 10 /zm R — 0.28). Tabulated values of 
R (A ) are given in Ref. 37 For sufficiently thin samples 
(ad < l) the transmittance measured as a function of the 
thickness of the sample can also be used to determine R (A ). 
Measurements carried out in this study showed that for elec- 
tron concentrations exceeding I0 18 cm' J and for A > 10 fim 
the reflectance decreases with increasing electron concentra- 
tion (see also Refs. 12 'and 15). However, for lower electron 
concentrations, R was found to be essentially independent of 
n . 

The transmission measurements, intended to provide a 
basis for determining the absorption coefficient, were car- 
ried out on thick samples and, with the exception of sample 
7, the condition ad > 1 was satisfied for wavelengths in the 
vicinity of 10 fim. 

The experimental results on the absorption coefficient 
for A= 10 /mi are given in Fig. 5 together with those.avail- 
able in the literature. For comparison, the theoretical values 



FIG 5 Theoretical and experimental values of the absorption coefficient as 
a function of electron concentration in GaAs at room temperature 


of the tolafabsorplion coefficient a t and the absorption com- 
ponents a opr a tmpl and a at . are also given m Fig. 5. As it lias 
been discussed above, these absorption coefficients were 
computed assuming no compensation. Thus, it is understan- 
dable that experimental values of a are greater than th** theo- 
retical ones Some of the experimental values of Ref ! i are 
smaller than those computed theoretically. Although the 
reason for this discrepancy is not dear, it should be noted 
that in the present study absorption coefficients smaller than 
those obtained theoretically were observed in samples exhib- 
iting signigicant carrier-concentration inhomogeneities. 

B. Characteristics of tree-carrier absorption 

As pointed out above, a comparison between the experi- 
mental and theoretical dependence of free-carrier absorp- 
tion on the electron concentration for samples of different 
degrees of compensation can be achieved by reducing the 
experimental absorption coefficient to zero com pen sa- 
tion[Eq. (13b)]. Such a reduction of the expenmen ta I results 
of Fig. 5 was carried out using the values of compensation 
ratio-determined for each sample from electron mobility 
(Table III), and the computed values of a xmp as-given in Ta- 
ble IL The results of a* xp versus electron concentration are 
given in-Fig. 6 together with the theoretical dependence of 
a t (n). The obtained agreement is very good, and it is worth 
emphasizing that no adjustable parameters were used in the 
calculation of the mobility, which provides the basis for the 
determination of 0 and of the free-carrier absorption 
coefficient. 

' The values of the compensation ratio determined from 
electron mobilities (Table III) have also been used to com- 
pare the experimental and theoretical dependence of the ab- 
sorption coefficient on the wavelength A. Typical results are 
given in Fig. 7. It is seen that very good agreement is ob- 



FIG. 6, Theoretical (.solid line) and experimental values ofihe absorption 
coelficient reduced to zero compensation (see text) as a function of electron 
concentration ui GuAsai room temperature; wavelength of radiation 2 = 10 
/an. 
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FIG. 7. The absorption coefficient as a function of wavelength Solid and 
dotted lines correspond to theoretical values, neglecting compensation, and 
dotted lines correspond to theoretical values obtained by taking into ac- 
count compensation as determined from electron mobility; •, ESI and A are 
experimental values obtained with samples 6, 3, and 2, respectively, of 
GaAs as in Table III. 


tained when the effect of compensation is quantitatively tak- 
en into account. It is also seen that for the lower electron 
concentration (n = 1.1 X 10 17 cm' 3 ) the compensation affects 
noticeably both the absolute value of a and its wavelength 
dependence (i.e., the slope d In a/d liU-). This behavior is not 
surprising since with increasing compensation the dominant 
role of.polar scattering is superseded by ionized-impurity 
scattering. These two scattering modes lead to different 
wavelength dependences of the absorption coefficient. 


C. Compensation ratio 


The present results have shown that the measured val- 
ues of electron mobility (Table lit) and the experimental 
characteristics of the free-carrier absorption (concentration 
and wavelength dependences) are satisfactorily explained by 
theory taking into account compensation. 

As discussed above, the value of the compensation ratio 
can be determined independently from the electron mobility 
and the electron concentration (utilizing Table II) and from 
the free-carrier absorption coefficient [utilizing E'q. (12b) 
and the values of absorption components given in Table II]. 
The values of 6 obtained by these methods are in very good 


agreement as shown in Tableau. Tne differences in values o 
6 are within expci imcntal error which is primarily due to 
mhompaeneitv of the materiar.Tn the case of intentionally 
compensated epitaxial layers (samples 4 and '■ in table III) 
the values of d determined (O.o7 and 0 / 6 , respectively) her 


are in good agreement with the nominal compensation ratio 


estimated from the concentration of impurities introduced 
during the growth (i e , 0.7 and 0.75, respectively). Thus, ii 
must be concluded that both procedures for determining the 
compensation ratio iead to reliable and consistent results. 

As pointed out above, the upper limit to which the pre- 
sent approach is applicable corresponds to electron (and im- 
purity) concentrations of about 3 X 10 ts cm" 5 . The lower con- 
centration limit is determined by the sensitivity of the room- 
temperature electron mobility and of the free-carrier.absorp- v 
tion to ionized lmpurites. These limits are apparent from 
Table II (and/or Figs. 4and 5). Thus, for a reliable determi- 
nation of the compensation ratio from electron mobility, the 
electron concentration should be greater than 10 !S cm • and 
in the case of free-carrier absorption the electron concentra- 
tion should be greater than 10 16 cm/ It should also be em- 
phasized that, for the reliable application of the present pro- 
cedures, GaAs should be reasonably homogeneous. Thus, it 
is desirable to supplement-measurements of mobility and/ or 
free-carrier absorption with measurements of the carrier 
concentration distribution on a microscale. Employing scan- 
ning ir absorption 35 for this.purpose, it also becomes possible 
(on the basis of the present study) to obtain a compensation 
ratio profile with a spatial resolution of about 20 /im. 

IV. SUMMARY 

Theoretical and experimental studies of room-tempera- 
ture electron mobility, and the free-carrier absorption coeffi- • 
cient have been earned out on n-type GaAs. It was, shown 
that for a total concentration of ionized impurities smaller 
than 10 15 cm" 3 , the room-temperature electron mobility is 
dominated by polar scattering (with about 10% contribution 
from acoustic-phonon and piezoelectric scattering) and it 
approaches a value of about 8000 cmW sec, in agreement 
with previous experimental and theoretical results. For im- 
purity concentrations exceeding 10 ,s cm' 3 , the mobility be- 
comes sensitive to ionized-impurity scattering, permitting 
the determination of the compensation ratio in the material. 

It was found, however, that the contribution from ionized 
impurity cannot be analytically separated from other scat- 
tering mechanisms since the application of Matthiessen’s 
rule. leads to significant errors. Thus, in determining, the 
compensation ratio, one must rely on the results of numeri- 
cal calculations. The calculated Values of electron mobility 
as a function of compensation ratio given in Table II can be 
readily used for the determination of the compensation ratio 
in w-type GaAs with electron concentrations ranging from 
10 15 to 3x 10 M cm" 3 . 

It was also shown that the free-carrier absorption coef- 
ficient is sensitive to compensation in GaAs for electron con- 
centrations exceeding 10 16 cm" 3 . An analysis of the present 
experimental results and those available in the literature 
showed that the characteristics of free-carrier absorption 
(concentration and wavelength dependences) can be satis- 
factorily explained only when compensation is quantitative- 
ly accounted for. 

A procedure for the determination of the compen ->:Uion 
ratio from the value of the free-carrier absorption coefficient 
for/t„= 10/rrn (or other wavelengths 8</U 12/z.m), and the 
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calculated values of the component absorption coefficients 
was Outlined For reasonably homogeneous materials, a 
good agreement was obtained between the values of compen- 
sation ratios obtained from electron mobility and from frec- 
Qarrier absorption it should' be pointed out that no titling 
parameters other than* the compensation ratio were used in 
comparing the experimental and theoretical results. 

It* was thus concluded that both the room-temperature 
electron mobility and the free-carrier absorption coefficient 
provide a reliable and consistent means for deternnningthe 
degree*of compensation (or the total concentration of ion- 
ized impurities) in n-type GaAs. In the form presented here, 
these methods can be used as practical procedures requiring 
relatively simple instrumentation. The electron concentra- 
tion and mobility can be determined by means of standard 
conductivity and Hall-effect measurements, taking the Hall 
factor r as unity; the error due to this assumption is estimat- 
ed to be about 10% for /i ~ 10 15 cm’ 3 and decreases with in- 
creasing carrier concentration. 

The free-carrier absorption method outlined here re- 
quires the determination of the absorption coefficient for one 
arbitrarily chosen wavelength in the spectral region 8<A< 12 
fim . Accordingly, relatively simple ir spectrometric instru- 
mentation is required (e.g., a C0 2 laser and suitable ir detec- 
tor). For sufficiently thick samples (satisfying the condition 
adj» 1) the error associated with minor changes in reflec- 
tance due to electron-concentration changes cambe-neglect- 
ed(for«<3x 10 ,s cm’ 3 ) and one can utilize the values ofR (>i) 
given in Ref. 36. When the electron concentration, the elec- 
tron mobility, and/or the free-carrier absorption coefficient 
are experimentally determined, the compensation ratio and 
the total number of ionized impurities can be determined 
from Table II and Eqs. (1 1) or (12). 
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APPENDIX 

. Equation (16) describing low-field electron, mobility 
has been obtained with the aid of the variational principle 
method 4 adapted to the case of the screened electron-opti- 
cal-phonon interaction/ 7 The change of the electron-distri- 
bution-function / can be written in the form 

(fT) = —Ak cost? ~-L (C), (Al) 

' \ at /coil ox 

where x = E/k 0 T y Eh the electron energy, /:cos<? is, the pro- 
jection of the electron- wave vector on the electric field direc- 
tion, and 

A=A(pc,T,m*) 


= 4.768 X 10 2s 




z t ( m*/m 0 ) l/ 2 /I \\ 

r 71 U. “ ej 


(A2) 


The operator L (C ) consists of four terms describing the dif- 
ferent scattering mechanisms considered at present 


L(C) = L t (C) 


C{£)( 1 


A \ * imp 


“h 


1 


4- 


4-) (A3) 

* pd / 


The relaxation times r for elastic scattering processes are 
given by the following: 

Screened ionized impurities, 32 


_L =2.415 W2 (xT)~ 3/2 

”imp 4 V m, / 


X 




Ax/a 


+4 x/a J 


(A4) 


Acoustic-phonon scattering through deformation 
potential? 3 


1 


=4.167 X 10 15 




pvj 


V 


m 0 


(A5) 


and piezoelectric (acoustic-mode) scattering, 34 

— = 1.052X 10 7 XA + i- )(- — )’ /2 . 

T pe , V c, C t J\x m„ J 


(A6) 


For screened optical-phonon (polar) scattering 27 
£ opl (C) = (fJf){N+ l)C.p. + a)S, - aR.T. - 4tt] 

- 2x{fJf){N + DOCS’. - aT.) + h(x- z,) 


X ( (NfJf)C-[(R. + a yS. - aR.T. - 4U.] 

- 2x'NC.(f./f)(S. r aT.)}, (A7) 

where •' 

f± —fo( x + z i)< C = = C(E + fico,), 

R ± = 2x-pa±z„ U+ — [*(x ±Zi)] ]/1 , 

*‘ ± -*o‘ ± ' p, c . 

0F POOR 

r*.+2u„\ ^Quality 


t _ = 


s.-J R * +2U * \ 

* \R„- 2 ’tr J 


and N= l/Qexp^^ — 1] isthe optical-phonon occupation 
number. Following the procedure developed in Ref. 4 one 
obtains the following expression for electron mobility: 


fi = 308.6[ ( V- — V /2 r 

LV ^0 A m o / 


1/2 


907 


J. Appl Phys., Vol 50, No. I February 1979 


Walukiewtcz etai. 


907 



37 - 



(A8) 


where the determinants have forms 


•P3/2, 3/2— 


0 

P$ n) 

p? /2) 

£o 3/2) 

do o 

d oi 

ft j 3/2) 

dio 

dn 


(A9-) 



(A10) 


The elements of the determinants are given by the integrals 


r E i/2 0 ^_ dE . • (All) 

• f Jo dE 

. (A12) 

, The functions <P r should represent a complete set of func- 
tions of electron.energy E. In the present computations, the 
set =- E —was chosen with values of r ands equal to 0, 1, * 

and 2, which at room temperature assure an accuracy of the 
mobility calculations within a few percent: 4 - 36 
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Concentration Distribution in Semiconductors 
by Scanning IR Absorption: Si 
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ABSTRACT 

A contactless method is reported for the quantitative measurement of 
spatial variations of the free carrier concentration m semiconductors with a 
resolution of about 20 and a sensitivity better than 1%. The method, is 
based on free carrier absorption in the infrared region where other absorption 
processes are negligible. In the experimental arrangement a CO 2 laser beam 
passes through a semiconductor wafer mounted on a stage with a calibrated 
x-y motion; 6 detection of the transmitted light through the desired area of 
the wafer about 20 ^m in diameter is achieved by employing an 1R micro- 
scope in conjunction with an IR detector. The principles and applications 
of the method are discussed using n- and p-type silicon. Carrier concentration 
profiles obtained by scanning IR absorption were tound to be in good agree- 
ment with those obtained by standard spreading resistance measurements. 


The determination of spatial variations of the carrier 
concentration in semiconductor materials continues 
to be of fundamental importance in device studies and 
applications and in studies of crystal growth and segre- 
gation. A number of methods have been developed for 
carrier concentration profiling utilizing the Schottky 
barrier (1), ESM (2), electrorefiectance (3), electro- 
luminescence (4), IR topography (5, 6), and spreading 
resistance (7). Most of these methods are qualitative in 
nature and their application is limited to specific con- 
centration ranges. The most commonly used method to 
obtain quantitative carrier distribution profiles is based 
on spreading resistance measurements. This method, 
however, is applicable only to elemental semiconduc- 
tors (pressure metal contacts with compound semicon- 
ductors are irreproducible) ; it introduces indentations 
on the surface and is very sensitive to surface prepara- 
tion conditions (7). Optical methods are most desir- 
able, as they require no electrical contacts (5,6). 
Thus far, striations in Ge (5), Si, and GaAs (6) have 
been observed through free-carrier IR absorption (8). 
Infrared topography is limited to a narrow carrier 
concentration range and is of a qualitative nature. 

In the present study IR absorption by free carriers 
in a scanning mode is employed for the quantitative 
determination of free-carrier concentration profiles in 
Si over a wide concentration range. 

Experimental Approach and Results 

Procedure : — Silicon samples of n- and p-type carrier 
concentrations ranging 10 15 -10 ia cm” 3 were cut from 
Czochralski-grown crystals. They were mechanically 
lapped to thicknesses ranging from 1 cm to 50 j±m. (de- 
pending on the carrier concentration) and Syton pol- 
ished. In order to assess the effect of surface prepa- 
ration on the reflection coefficient, final polishing was 
also carried out with 0.3 or 0.03 /un particle diameter 
abrasive. 

• Electrochemical Society Active Member. 
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Measurements of infrared transmission in the spec- 
tral region of 2.5-25 Am were performed in a Fourier 
spectrometer. The values of the absorption coefficient, 
«, were obtained from the measured transmittance T 
using the expression 

= J_ 

Jo 1 - R2 e -2a* 

where I is the intensity of the transmitted light, I 0 is 
the intensity of the incident light, R is the reflection co- 
efficient, and z is the thickness of the sample; the spec- 
tral dependence and values of R employed in this study 
were taken from Ref. (9). Equation [1] was also used 
in some instances to obtain the value of R by mea- 
suring the transmittance as a function of sample thick- 
ness (10). With this procedure the reflection coeffi- 
cient at 10 /im was obtained for sample with different 
carrier concentrations and different surface prepara- 
tions. 

Dependence of absorption coefficient on wavelength 
and carrier concentration . — The measured absorption 
coefficient of n- 'and p-type Si as^a* function of wave- 
length is presented in Fig. l~and £ for various concen- 
trations. It is seen that, in general, the absorption co- 
efficient increases with increasing wavelength; This in- 
crease is primarily due to free-carrier absorption (8) ; 
infrared radiation is absorbed by free carriers in the 
valence or in the conduction band when scattering 
mechanisms (by phonons or ionized impurities) pro- 
vide conservation of momentum. For small carrier 
concentrations, other absorption processes leading to 
strong absorption peaks are superimposed on free car- 
rier absorption. 

In the vicinity of 10 fim (outside of the, oxygen ab- 
sorption band) even for small-carrier concentrations, 
the free-carrier absorption dominates all other ab- 
sorption processes. The absorption coefficient of n- ana 
p-type Si taken from Fig. 1 and 2 for radiation with a 
wavelength of 10 is shown in Fig. 3 and 4, re- 
spectively, as a function of free-carrier concentration, 
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X [/tm] 

Fig, 1. Dependences? absorption coefficient, a , on wavelength of 
radiation for Si samples with different electron concentrations: 
1, U X 10^; 2, 7.6 X 10 13 ; 3, 5.1 X 10 16 ; 4, 2 X W 7 ; 5, 
7.2 X 1G 17 ; 6 , 2.9 X 1018; 7 r 3.6 X I0 18 ; and 8 , 7.1 X 10 18 
cm“ 3 . 



X f/imj 

Fig. 2. Dependence of absorption coefficient on wavelength of 
radiation, a, for Si samples with different hole concentrations: 
1, 7.4 X 10 15 ; 2, 3.3 X 10 16 ; 3, 9.1 X 10 18 ; 4, 1.7 X 10 17 ; 5, 

5.2 X 10 17 ; 6 , 7.4 X 10* 7 ; 7, 1.4 X 10 18 ; 8 , 5.0 X 10 18 ; 9, 

1.2 X IQ 19 ; and 10, 1.7 X 10 w cm"» " ^ 


together with the data reported in the literature (11- 
14). In n-^type Si, the absorption coefficient increases 
linearly with increasing electron concentration rang- 
ing from 10 15 to about 5 x IQ 17 cm -3 ; for higher elec- 
tron concentrations the absorption coefficient exhibits 
a slightly overlinear dependence on the carrier con- 
centration. This behavior is associated with an in- 
creased contribution of impurity scattering (8) to free- 
carrier absorption. 

In p-type Si (Fig. 4) the absorption coefficient ex- 
hibits a linear dependence on hole concentration ,in the 
range of 3 X 10i 6 -2 x 10 19 cm“ 3 . For smaller concen- 
trations the absorption coefficient decreases, but less 
than expected for free-carrier absorption. 

Dependence of transmittance on carrier concentra- 
tion.— The dependence of the' absorption coefficient on 
the carrier concentration (for >. = 10 ^m) can be em- 
ployed for the determination of the absolute* value of 
the carrier concentration and carrier concentration 
fluctuations. It is convenient to introduce a cross sec- 
tion of absorption, Kn (P , describing the sensitivity of 
the absorption coefficient changes to carrier concen- 
tration changes, as follows 

[ 2 ] 

n,p 

For n-type Si, K n = lx 10” 16 cm” 3 (for electron con- 
centrations up to 5 X 10 17 cm” 3 ) and increases to 2 x 
10“ 16 cm” 3 as the concentration increases to 2 x 10 19 
cm” 3 (Fig. 3). For p-type Si, — 2 x 10” 10 cm” 3 
for concentrations higher than. 3 X 10 IG cm 3 ; * or small 
chances in carrier concentration (An,Ap), Aa =K n . P An, 


a 



Fig. 3. Dependence of the absorption coefficient on electron 
concentration for radiation with 10 wavelength: e present 
work; A Ref. (11), Q 03). 
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Hole concentration, [cm' 3 ] 


Fig. 4. Dependence ohabsorption coefficient on hole concentration 
for radiation with 10 pm wavelength: # present work; A Ref. (14), 

O Ref. (12). 

Ap. Changes in concentration .should lead to relative 
changes in transmittance, and if changes in ft are 
negligible, then according to Eq. [1], for An,Ap « n,p, 
one obtains 

1 AT 

- “ V— = X An, Ap [3] 

z T av 

where T av is the average value of transmittance along 
the wafer. Thus, variations in carrier concentration 
can be determined from variations in transmittance, 
providing corresponding variations in R are insignifi- 
cant. 

Reflection coefficient— The reflection coefficient be- 
comes sensitive to carrier concentration variations only 
when the wavelength of the radiation is near the 
wavelength at which the minimum in reflection due 
to plasma absorption occurs (9). For Si with mod- 
erate carrier concentrations the minimum in reflec- 
tion occurs at far longer wavelengths than that of the 
radiation analyzed in this study (10 pm). For n-type 
silicon this minimum begins to shift appreciably at 
carrier concentrations above 10 18 cm" 3 , and reaches 
a wavelength of 22 ^m at 5 x 10 18 cm" 3 ; in p-type 
Si the minimum shifts to 16 ^m at a carrier concen- 
tration of 1.2 x 10 19 cm" 3 (15). 

Values of R were determined as a function of car- 
rier concentration employing Eq. [1], as outlined 
above. It was found that R remained constant (0.29) 
for concentrations up to 1 x 10 18 cm" 3 in the case of 
n-type Si and for concentrations up to 2 x 10 18 cm -3 
in the case of p-type Si; for higher concentrations ft 
decreased to 0.26 in going to n = 5 x 10 18 cm" 3 and 
to 0.19 for p = 1.2 x 10 19 cm" 3 . However, these 
changes of ft were found to be orders of magnitude 
-smaller than for corresponding changes in the absorp- 
tion coefficient Thus, in going from p = 0.9 x 10 19 
cnr* 3 to p = 1 X 10 19 cm" 3 , the contribution of the 
change of' ft to the change of transmittance was found 
to be an order of magnitude smaller than the contribu- 
tion of the absorption coefficient for samples where 
az = 2. Accordingly, for samples where oz > 2, 'varia- 
tions of ft are not taken into account in analyzing 
variations in transmittance. 


The effects of variations in surface preparation on ft 
were assessed by determining ft for samples having the 
same carrier concentration but having been polished 
with an 0.3 and 0.03 pm particle diam abrasive and by 
Syton. The value of ft was found to be the same, 
within experimental error (about 5%), for samples 
polished with 0.03 p m particle abrasive and for those 
polished with Syton. For samples polished with 0.3 
^m particle diam abrasive the value of ft was about 
10% greater than that for samples polished with Syton. 
This variation of ft for samples where az r: 2 con- 
tributes about 3% error to the values obtained for 
and for carrier concentration. Accordingly, variations 
in ft due to surface preparation variations need not 
be taken into account in analyzing variations in trans- 
mittance. 

Scanning 1R absorption measurements . — Scanning 
IK absorption measurements were carried out employ- 
ing an experimental arrangement shown schematically 
in Fig, 5. A parallel beam 1 cm in diam from a CO> 
laser (tunable between 9.16 and 11.02 pm) is trans- 
mitted through the sample positioned on a 'stage with 
x-y motion and scanning rate ranging from 50 to 10 3 
/jtm/min. The light transmitted through the desired 
area of the sample is focused on a detector (bolometer) 
through the optics of an IR microscope, and the signal 
from the detector, which is proportional to the inten- 
sity of radiation, is amplified and' recorded. Using a 
74 x objective, the image of an area of about 20 pm 
diam is- obtained. 

In the present study a radiation power of about 0.2W 
was used to prevent heating of the sample. The aver- 
age transmittance (T av ) was determined from the re- 
corded intensity of the laser beam with and without 
the sample on the microscope stage. Changes in trans- 
mittance (AT) were recorded as the stage moved at 
a constant rate in the x or y direction. In the pres- 
ent arrangement changes in transmittance A T/T aV i 
with an accuracy of about 0.1%, and a spatial resolu- 
tion of about 20 p m were determined. On. the basis of 
Eq. [3] (taking the value of K n , P from Fig. 3 or 4), 
transmittance changes of 0.1%, correspond to fluctua- 
tions in carrier concentration of An = 10 l3 /z cm -3 
and Ap 5 x 10 13 /z cm" 3 ; thus, for a typical wafer 
with a thickness of 200 pi n, An = 5 X 10 14 and Ap = 
2.5 x 10“ cm" 3 

A typical example of carrier concentration fluctua- 
tions obtained by scanning IR absorption for an n-type 
Si sample, 100 /xm thick, with an average carrier con- 
centration of 3 x 10 18 cm" 3 (cut from Czochralski- 
grown crystal along the growth direction) is shown in 
Fig. 6. The observed carrier fluctuations of about 15% 
are typical for the commercially grown single crystals 
by the Czochralski method; these carrier concentra- 
tion fluctuations reflect turbulent convection-controlled 
segregation. 

Comparison of scanning IR absorption with spread- 
ing resistance measurements.— Scanning IR absorption 
was directly compared with spreading resistance mea- 
surements by obtaining carrier concentration pro- 
files on the same Si sample (average electron concen- 
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Fig. 5. Experimental arrangements for scanning absorption 
measurements. 
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Fig. 6. Carrier concentration fluctuation along the growth direc- 
tion in Czochralski-grown n-type Si. 

tration 2 X 10 18 cm” 3 ) employing both methods. To 
carry out spreading resistance measurements (one- 
point probe) a uniform ohmic back contact was made 
-to the sample with an indium layer and a scan was 
performed at 10 intervals; the ohmic contact was 
lapped off and the sample polished; subsequently, scan- 
ning IR absorption measurements were carried out 50 
/*m from either side of the spreading resistance scan. 
The results are shown in Fig. 7. It is seen that the 
results obtained by the two methods are in very good 
qualitative agreement. The observed differences such 
as changes in the shape and position of carrier concen- 
tration maxima are attributed to the fact that the scans 
are 50 /on apart and, thus, reflect radial variations in 
carrier concentration. In addition, the volume probed by 
the two methods is somewhat different; thus, in spread- 



O 500 JOOO 1500 

DISTANCE [>mj 


Fig. 7. Comparison of carrier concentration fluctuations in re- 
type Czochraiski-grown Si (average concentration 3 X 10 1S cm -3 ) 
by spreading^ resistance measurements (SR) and by scanning ab- 
sorption, 50 /xm from either side of the spreading resistance scan. 


ing resistance measurements the volume probed is 
semispherical with a diameter of about 20 Jim in the 
immediate vicinity of the pressure contact, whereas m 
JR absorption the volume probed is cylindrical with a 
diameter of about 15 and a length equal to the sam- 
ple thickness. 

Summary and Conclusions 

Infrared absorption by free carriers was employed 
for the quantitative determination of the carrier con- 
centration on the microscale m n and p Si. Employing 
radiation with a wavelength of 10 /tm, the method was 
found to be applicable to concentrations ranging from 
10 16 to 10 I() cm” 3 with a sensitivity of about 5 x 10 w 
cm” 3 and a spatial resolution of about 15 ^m. 

Carrier concentration profiles obtained with this 
method and with spreading resistance measurements 
were found to be in very good agreement. Thus, scan- 
ning IR absorption provides a fast, sensitive, contact- 
less, and reliable method for the determination of 
carrier concentration on a microscale in Si over a wide 
range of carrier concentrations. The application of this 
method for obtaining carrier concentration profiles in 
compound semiconductors is currently being, studied. 
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On the basis of mass-transport principles, a theoretical model of elect roepit axial growth — current - 
controlled liquid-phase epitaxy— was developed which defines the contribution of the Peltier effect (at the 
sohd-solution interface) and thatj)f solute elect romigration to the overall growth process According. to the* 
model, the contribution of elect romigration to growth is dominant in the absence of convection in the 
solution, whereas the contribution of the Peltier effect can be dominant in the presence of convection On 
the basis of the model, expressions were derived which relate quantitatively the growth velocity to growth 
parameters. The model was found to be in excellent agreement with extensive experimental data on the 
electroepitaxial growth of GaAs from a Ga«As solution 


PACS numbers- 68.55 -t-b, 81.I0.Dn, 68 45 -v, 61 50 Cj 


INTRODUCTION 

The electronic-device structures of GaAs and of all oth- 
er III-V compounds are based totally on epitaxial layers 
grownon Bulk single-crystal substrates. For most types of 
devices, epitaxial layers are presently grown from liquid so- 
lutions. In liquid-phase epitaxy (LPE) the substrate (e.g., 
GaAs)„is thermally equilibrated with the solution [e.g., 
Ga(rich)-As solution] and growth takes place upon supersa- 
turation of the solution brought about by a temperature 
decrease.. 

While LPE has been highly developed.in recent years, 
precise control of the microscopic growth velocity (and thus 
control of segregation and defect structure) is still a problem. 
It has been recognized for many years that control of the 
growth velocity and of segregation could be achieved 
through passage of electric current across the growth inter- 
face since it induces Peltier cooling (or heating)in the imme- 
diate vicinity of the interface 1 - 2 and also causes electromigra- 
tion of the solute to the growth interface. 3 " 8 Although 
approaches to the growth of bulk single crystal from the melt 
have been proposed, 1 - 2 - 4 they have not been reduced to prac- 
tice because of convective instabilities in the melt and other 
complexities introduced ;by Joule heating. 7 

Current-controlled grpwth has been successfully car- 
ried out in recent years 9 in an LPE configuration where Joule 
heating presents no significant difficulties because the di- 
mensions of the substrate and the solution are relatively 
small. Electroepitaxial layers were first obtained using the 
InSb system 9 ; growth was initiated and sustained by passing 
electric current through the substrate-solutiominterface 
while the temperature of the overall growth system was 
maintained constant. Electroepitaxydias since been em- 
ployed for jhe growth of GaAs, ltM7 InP, ,s GaAIAs, 19 " 23 and 
garnet layers. 24 

Electroepitaxy is still in the early stages of develop- 
ment. However, its advantages m achieving controlled dop- 
ing; 10 * 1 113 ' 14 improved surface. morphology, and defect struc- 


J> On leave from the Institute of Physics, Polish Academy of Sciences, Wai* 
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ture, 16,17 improved electronic structures, 20 13 and.in studying 
growth and segregation 14 have already been demonstrated. 

Attempts to explain the electroepitaxial process have 
been, thus far, qualitative in nature and limited in scope. It 
was originally proposed 9 that .Peltier cooling at the growth 
interface was responsible for supersaturation and electroepi- 
taxialgrowth. Subsequently, it was demonstrated that under 
certain conditions electromigration constitutes the domi- 
nant contribution to supersaturation, 12 while in other studies 
it was shown that Peltier cooling dominates electroepitaxy. 21 

In the present study a comprehensive theoretical model 
of electroepitaxy is developed which relates; on a quantita- 
tive basis, the contribution of Peltier cooling and electromi- 
gration to the growth kinetics; the model also relates quanti- 
tatively the growth velocity to growth parameters and it is 
successfully used for the analysis of extensive experimental 
results. In part II of this investigation a theoretical model of 
segregation in electroepitaxy will be developed. 

PROCESS AND PHENOMENA 
The process 

Electroepitaxy utilizes the standard LPE cpnfiguration 
modified to permit passage of electric current through the 
solutiomsubstrate interface: 15 A schematic representation of 
a typical growth cell is shown schematically in Fig. 1. After 



TEMPERATURE - 



FIG. 1. Schematic representation of growth cell used for electroepitaxial 
growth of GaAs from Ga-As solution and of the temperature profile result- 
ing from Peltier cooling at t he susbt rate-solution interface and from Peltier 
heating at the substrate-clecti ic contact intei face. 
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(a) (b) (c) 


PELTIER COOLING ELECTRO MIGRATION PELTIER COOLING PLUS 



In Ga- As solutions, eleclromigration of As-spectes is 
anode directed; thus, when the substrate has a positive polar- 
ity, the solution becomes supersaturated with As at the sub- 
strate-solution interface leading to epitaxial growth. 25 As 
shown schematically in Fig. 2(b). eleclromigration causes 
(in a thermally equilibrated solution) transition /f 0 — >A l , he., 
supersaturation, which in turn leads to, transition A X -+A 0i 
i.e., precipitation (solidification) of GaAs; in T,C coordi- 
nates, these transitions can be expressed as 
(TSCoyAT^Co+A C l )->(r 0 ,C 0 ). 


FIG 2. Schematic representation of the pertinent part of the phase diagram 
of an A m -B v system and transition involved in electroepitaxial growth. 


thermal equilibration,. growth is initiated and sustained by 
passing electric current through the substrate-solution inter- 
face in the appropriate direction (see below), while the over- 
all temperature of the system is maintained constant. 


Peltier cooling 


The substrate and solution, being different electrical 
conductors, have different thermoelectric coefficients. Thus, 
flow of electric current across their interface is accompanied 
by absorption of heat or evolution of heat, depending on the 
current direction. The magnitude of this heat, Q, is propor- 
tional to the difference in P.eltier coefficients, ir pi and the 
current density For A lll B y semiconductor com- 

pounds and the temperature range used in electroepitaxy, Q 
is of the order of 1 W/cm 2 at a current density of about 10 
A/cm 2 : In an equilibrated isothermal LPE system the ab- 
sorption of this heat decreases the interface temperature by 
AT py which is typically of the order' of O.Lor 1 °C and, thus, 
induces supersaturation leading to epitaxial growth.. 

. Using the pertinent part of a schematic-phase diagram 
of the compound A lll B v [Fig. 2(a)], Peltier-induced LPE can 
be illustrated by. the vertical transition^— super- 
. cooling; and the subsequent horizontal transition^ 1 — >A lt 
which leads to precipitation (solidification of compound 
A m B v ). Using temperature-composition coordinates (T,C), 
the process can-be expressed as (TLC*)— ^(r 0 — AT p ,C 0 ) 
-+(T 0 —AT p ,Ci). Accordingly, in the following treatment of 
^growth induced by Peltier cooling, the supersaturation 
' Cj — C 0 wilf'be approximated as ' - ' 


C,-Co= 



T 0 - 

y 

L 


where L designates liquidus. 


( 1 ) 


Eleclromigration 

The solutions^ elements from groups III-V used in 
electroepitaxy are metallic conductors; essentially they ex- 
hibit no ionic contribution to the current. However, in these 
solutions, as in other liquid metals, electromigration takes 
place due to electron-momentum exchange and/or electro- 
static field-forces. 8 Under the electric field E induced by the 
currenFflow, species in solution migrate with a velocity 
v=fiE, where jx is the mobility. 


Simultaneous Peltier cooling and 
electromigration 

When both Peltier cooling and electromigration are sig- 
nificant, the following transitions are involved in.the growth 
process [Fig. 2(c)]: the diagonal supersaturation transition 
A 0 —*A\ i.e., ( T 0t Co)^(T 0 -A T p , C 0 +A C) and the horizon- 
tal precipitation (solidification) transition A '—+A lt i.e., 

( T 0 -A T p ,Co+AC)^(T or A T py C t ] ). 

In treating the solidification (precipitation) transitions 
as strictly' isothermal, it must be assumed that the contribu- 
tion of the heat of solidification is negligible. This assump- 
tion is justified for the commonly encountered rates of solidi- 
fication in LPE which are_of the order of 1 /zm/min. 


MODEL OF ELECTROEPITAXY 
Growth kinetics ... 

In the model of electroepitaxy presented in this paper it 
is assumed, as in LPE, 26 ' 27 that the solute transported to the 
advancing growth interface is removed from the solution 
only through epitaxial growth on the substrate. Transport of 
solute due to temperature gradients in the solution is neglect- 
ed but solute transport due to convective flow in the solution 
is taken into consideration/In certain approximations and in 
calculations used in the present treatment, electroepitaxy of 
GaAs from Ga-As solutions is used as a specific-case. 

The widely adopted isothermal diffusion treatment of 
LPE 26 27 is here extended to include solute transport by elec- 
tromigration which is controlled by the mobility /z of the 
migrating solute and the electric field is in the solution: thus, 


n d 2 C dc „dC dC 

dx 2 dx dx dt 


( 2 ) 


where D is the diffusion coefficient of the solute in the solu- 
tion, C is the solute concentration, x is the distance from the 
advancing growth interface, u is the growth velocity, and t is 
the time. (The sign of the electromigration term fiE dC/dx is 
determined by the direction of the current flow.) 

The following boundary conditions are considered ap- 
plicable in the present case (the notation is as in Fig. 2; the 
concentration in the solid is € s and the thickness of the solute 
boundary layer is 5). 


D 


d£ 

dx 


—fiECi = v(Ci — 
o 


C s ), 


(2a) 
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at 0, 

c=c„ 

for all .v, (2b) 

at t > 0 t 

c=c. 

for* = co (absence of convection), (2c) 

or 



at f>0; 

C=C, 

for * > <5 (presence of convection), (2d) 

at ^>0 and x~ 0, 

C—C { (growth follows the liquidus line). 


(2e) 


In solving the transport equation (2), the term 
v(dC/dx) is neglected as- being very small compared with the 
diffusion and the electromigration terms; for example, in 
electroepitaxy of GaAs, fiE is usually one order of magni- 
tude greater Than v; in addition, it is assumed that the con- 
centration of the solute for either *= oo or x— 8 remains 
constant during growth; this assumption is considered valid 
when feed material (on top of the solution) is present or when 
the thickness of the grown layer is very small; in the latter 
case.the maximum time of growth is.limited and in common 
LPE is of the order of h 2 /D } where h is the height of the 
solution. 28 

The procedure for solving Eq. (2) is presented in Ap- 
pendix A. The expressions obtained for the growth velocity 
of electroepitaxy are given below: 

In the absence of convection (5= oo): 

p AT p dC 1 (D\n 

■ C s -C t dT\ L \irt ') 

exp[-(p.Et) 2 -MDt] _ C, fTl 

sdc[-/nEt/2(Dt) l/1 ] C s —Ci ’ 

in the presence of convection (finite boundary layer 
thickness) 

AT p dC I fDyn 

C‘ M -Ci dT\i\irtJ 

exp[—(/iEt) 2 /4Pt ] 

erfc[ — fiEt/ 2{Pt ) !/2 ] — erfc [(S—^iEt )/2,(Dt ) ,/2 ] 


— fiE 


C, 

C s — Ci ' 


(4) 


where A T p is the temperature decrease at the interface due to 
Peltier cooling. 

Equations (3) and (4) can be rewritten in a simplified 
general form where the contribution to the growth velocity 
by the temperature decrease at the interface (Peltier cooling) 
and.by electromigration are designated as u r and v E , 
respectively: 

u=v T f k (E,5,t)+v £ , (5) 

where 


AT p dC_ 1 /AV /2 

C.-C.^Llwr/ 


(6a) 


v E —fiE 



(6b) 


in the absence of convection: 

exp,[ -(/iEt f/ADt 


A erfc[ —fiEt/ 2(Pt) l/1 ] 

and in the presence of.convection 


f k =m,8,t) 


<^'®00R «° P 


(6c) 

<1$ 


exp.[ — (uE/ ) 2 /4Zh ] 


’ erfc[ —/iEt/ 2(i)/ ) ,/2 3 -erfcl )/2(Dt)' 


(6d) 

Note that for£— >-oo,/j(£,5,r )—>f x (E,t ). 

The asymptotic’behavior of/isas follows: 

f,(E,t)^(7rt/py /2 fiE for y.E- (t ) ,n 42(P ) ! ' : , (7a) 

for fiE (t 2(D ) 1/2 , (7b) 

f 2 (JE,t,8)~(irDt y n /8 for 8 l /AD4t4,AD/p 2 E l . (7c) 

It is of interest to point out that in view of the above 
asymptotic behavior of f, Eqs. (3) and (4) derived here, re- 
duce for E — 0 and A T p — A T to the standard expressions 
for diffusion-controlled growth 26 ”: 


AT dC 1 

C s -C, dTlAirt J ’ 

(8a) 

AT dC\ D_ 
1 C s —Ci dT \l8 

(8b) 


Absence of convection 


In electroepitaxy from nonmetallic solutions — e.g., 
electroepitaxial growth of garnets 24 — the electric field can be 
of the order of V/cm so that after a time period of a few 
seconds the inequality fiE-(t ) ,/2 >2(Z> ) 1/2 is satisfied. Accord- 
ingly-, approximation (7a) is valid and the time dependence 
of v in Eq. (5) cancels out through v r [Eq. (6b)] and /(£,?); 
thus, Eq. (3) reduces to 


fiE (AT p dC 
C s -Q \ 2 dT 


•+C t 

L 


)■ 


(9) 


In such a case, electroepitaxy proceeds at a constant growth 
velocity for a constant current density; actually, the growth 
velocity can be proportional to J 2 , since both E and A T p are 
proportional to /; diffusion plays no role in the growth 
process. 

Oiuthe other hand, in metallic-type Ga-As solutions a 
typical value of pE for 10 A/cm 2 is* of the order of 10" 5 
cm/sec 1229 and the diffusion constant of As (at 800 °C) is 
about 6x 10“ 5 cm 2 /sec J0 ; thus, for time periods less than 
about 10 h (/< 2X 10 6 sec) approximation (7b) is valid. Ac- 
cordingly, Eq. (5) reduces to 


v=v r + v 

C s -Q dTUKirt J C s -C x 


( 10 ) 
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FIG. 3. Growth (or dissolution) velocity transients of GaAs calculated 
from Eq. (10) from a Ga-As solution taking r o =800 *C t D— 6x 10* J 
cm 2 /sec,/i£= 10' 5 cm/sec, and A T p = — 3 °C for Peltier cooling and 
AT p — +3 *C for Peltier heating at the substrate-solution interface The 
contribution of the Peltier effect and of electromigration to the overall 
growth (or dissolution) are indicated. 


The first term in Eq. (10) represents the contribution of the 
Peltier effect, whereas the second term represents the contri- 
bution of electromigration to the growth velocity; both 
terms depend linearly on the current density, through A T p 
and E t respectively. The relative contribution of these terms 
to electroepitaxy of GaAs from a Ga-As solution at 800 °C is 
shown in Fig. 3; a A T p value of — 3 °C is assumed and fiE is 
taken equal to 10“ 5 cm/sec for a current density of 10 A/cm 2 . 
It is seen that after 1 h of growth the contribution of Peltier 
cooling to the growth rate reduces to about 15%; it is also 
seen that even in the presence of Peltier heating at the inter- 
face, electroepitaxial growth does take place, i.e., the contri- 
bution of electromigration to growth dominates over disso- 
lution of the substrate due to Peltier heating. 


On the basis of the present model the salient features of 
electroepitaxial growth in the absence of convection can be 
summarized as follows: (a) In an isothermal substrate-solu- 
tion system electroepitaxial growth can be initiated and sus- 
tained through solute electromigration to the interface; the 
contribution of Peltier cooling at the interface to growth is 
relatively small, (b) In metallic-type solutions, as in those of 
A ni B v compounds, the growth velocity is proportional to 
the current density since the electric field in the solution and 
the Peltier cooling are proportional to the current density 
(Eq. 10). (c) The contribution of the Peltier effectto the 
growth velocity decreases with time as \/{t ) !/2 , which is a 
typical behavior of diffusion-controlled processes; the con- 
tribution of solute electromigration to the growth velocity is 
time independent, (d) The temperature dependence of the 
growth velocity is primarily dictated by the phase diagram of 
the system (e.g.. Fig. 2). The Pel tier-effect contribution is 
also enhanced to some extent by temperature through D. 
The dependence of ir p and/t on temperature is relatively 
weak. 


Presence of convection 

Convection in the solution leads to the formation of a 
solute boundary layer of thickness <5. Approximation l ~c) is 
readili satisfied after a few minutes of grow ^accordingly, 
Eq. (3) reduces to 


C t — Ci dT' l8 c,-c, 


( 11 ) 


The presence of a boundary layer significantly increases 
the contribution of the Pei iter effect to electroepitaxial 
growth, as illustrated in Fig. 4(a), where the growth velocity. 



FIG. 4(d). Grow th % elocity of GaAs from a Ga-As solution calculated from 
Eq. (1 i) for various values of the solute-boundary layer S, taking 
T 0 — 800 "C. D=6 X 10 ' cmVscc, and //£ = 10 ' cm/sec. Peltier cooling re- 
sults when an /Mype substrate has positive, polarity: v h is the contribution 
ofclectromigration to growth velocity, (b) Growth velocity (solid curves) or 
dissolution n elocity (dashed curves) of GaAs calculated from Eq f 11) for 
different \ alues of boundary-layer thickness S. as a function of Peltier heat- 
ing (p-tv pc substrate with positive polarity) taking T c , D t a»d//£asin ta). iy 
is the contribution of electromigration to the growth velocity. 
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as calculated from Ecj. (1 1), is plotted as a function of —AT p 
(Peltier cooling) for various values of <5. The parameters used 
in these calculations arc the same as those used for Pig. 3, i.c., 
they correspond to GaAs growth from a Ga-As solution at 
800 *C under 10 A/em\ The contribution of electromigra- 
tion (which is independent of A T p ) to the growth velocity is 
about 0.3/*m/mm (as in Fig. 3). For a given (small) \alue of 
5 the growth velocity increases rapidly with increasing 
Peltier cooling <J T p ) at the interface, 'and for a given 
amount of Peltier cooling the rate increases rapidly with de- 
creasing 5 (increasing. convective flow) However, for large 
values of $ ( > 10 H cm) and A T p values of the order of 
—0.1 °C, the growth rate is controlled by electromigration. 

In Fig. 4(b) the growth (or dissolution) velocity calcu- 
lated from Eq. ( 1 1) [as in Fig. 4(a)] is plotted as a function of 
Peltier heating for various values of 5. It is seen that growth 
does take place for relatively small amounts of Peltier heat- 
ing and/or large values of <5. Dissolution of the substrate is 
pronounced for small values of 5 and particularly as, A T p 
exceeds about 0.15 °C. 

Thus, in the presence of convection the salient features 
of electroepitaxial growth are the following: (a) The contri- 
bution of the Peltier effect can dominate that of electromi- 
gration, particularly for small values of S and larger values of 
T p . (b) After a period of growth r>5 2 /4D, the growth is 
independent of time, (c) The dependence of the growth ve- 
locity on' temperature is as outlined in the absence of 
convection. 

Maximum thickness of epitaxial layer 

When no feed material is used in the growth system, 
electroepitaxial growth proceeds up to a certain thickness 
and then stops due to the gradual depletion of the solute in 
the solution. The.maximum thickness of the epitaxial layers, 
Z max < is related to fundamental growth parameters; since it 
can be readily measured experimentally, it contitutes andm- 
portant means for the characterizatiomof the growth 
process. 

The amount of solute incorporated in the epitaxial layer 
represents the decrease of solute concentration in the solu- 
tion;^ is, thus, apparent that 

z = Coh ~ £ C (*) dx, (.12) 

where x—0 at the growing interface, h is the height of the 
solution, and C (,r) is the solute concentration in the solution 
when the growth process stops; it is assumed that Z nu ^<h . 

Electroepitaxial growth stops when the net solute flow 
in the solution equals zero; thus, Eq. (1) reduces to a steady- 
state expression 

D „ ;i £C= 0. (13) 

dx 

the boundary conditions are 


C ll ^C l -(fiE/D)SC l . 


Thus, the solution of Eq. (13) becomes 


CW-C,*(d-C rf )X 


zxp(/iEx/D)— 1 
1 — exp (f/ES/D) 


(14) 


for .\<£ ' (15) 

and C (x)—C d , for a ><5 

For GaAs, fijtE/D^ 1. for h < 1 cm. Thus, for solution 
heights not significantly more than 1 cm, Eq. (15) reduces to 

C (jc) - C, = ( C d — C,)(x/8) for x<S ( 1 6) 

and 

C ( s ) =C d , for %><S. 

From Eqs./12), (14), and (16), Z max is readily obtained: 
thus, for <5= A (absence of convection) 

<i,a) 

and for <5</z (presence of convection) 

Z max ~ C °~ Cl h + f^£ls(h- — ). (17b) 

C s _ DC S \ 2/ 

It is apparent that in Eqs. (17a) and (17b) the first term 
represents the contribution of the Peltier effect and the sec- 
ond term the contribution of electromigration. 

In Fig. 5 the Z max values obtained from Eq. (17b) are 
plotted as a function of —A T p for various values of 5. The 
calculations were carried out for GaAs growth from a Ga-As 
solution at 800 °C, for A = 1 cm ,fiE= 10' 5 cm/sec, and 
D=6x 10" s cmVsec (C 0 and Q values are taken from the 
GaAs phase diagram). 32 Here again, for small values of 



C—C u for_v»0 (at rhe interface), 

C— C lt y for x>S* 

The value of C (i is obtained from Eq. (11) taking 0: 


FIG. 5. Maximum layer tluci ness of GaAs grown from a Ga-*V> solution 
calculated from Eq ( Vb) as a 1 unci ion of Pel her cooling at the ero a tii 
mtcrface for various values or solulc-ooimdary -layer lliickn^s «), i /I 

D y and uE as in Fig 4(a) 
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convection), the contribution of the Peltier effect to the 
maximum thickness of the epitaxial layer is verv small- the 
contribution of the Peltier effect becomes dominant in the 
presence of convection (small value of 5). 


Criteria for the contribution of Peltier cooling 
and of electromigraiion: critical Peltier cooling 
and critical substrate thickness 

According to the model developed and discussed above, 
there are two effects operating in electroepitaxy, Peltier cool- 
ing at the growth interface and electromigration of solute 
species to the interface (or solution species awav from the 
interface). The relative contribution of these effects depends 
on the experimental conditions. 


On the basis of the present model, quantitative criteria 
can be established for the determination of the relative con- 
tribution of the Peltier effect and electromigration to 
electroepitaxy. . • . 


It is convenient to introduce a critical value of tempera- 
ture decrease at the interface, A T c , for which the contribu- 
tion of the Peltier effect and electromigration to the growth 
velocity are equal; i.e., considering Eq. (5) 


\ v rf k \ = \ v E \. 

On the basis of Eqs. (6a) and (6b) 
AT - A-gCi dT j /£\w 

" f k dclAnJ 


(IS) 

(19a) 


In the presence of convection and fo'r t>8VD, Eq. Vl9a) re- 
duces to 


A T' -fxECt —I 
dC I 


8 _ 

lD 


(19b) 


Thus, for small amounts of Peltier cooling, A T p ^A T c , elec- 
tromigration dominates, whereas, for large amounts of cool- 
ing, A T p ^>A T the Peltier effect dominates the growth 
process. 


d — 2 l l P h ' C < clT I 

Vifx dC\ L \DJ 

In the presence of convection Eq. (21a) reduces to 

d ~mpkC, dT 5 

tt p dC lD 


(21a) 


(21b) 


where for d>d c the Peltier effect dominates and for d <d 
electromigration dominates the growth process. 

Thus, with thick substrates and in the presence ofsig- 
mficant convective flow, electroepitaxial growth is dominat- 
ed by Peltier cooling; on the other hand, with thin substrates 
and in the absence of significant convective flow, elect romi- 
gration dominates electroepitaxial growth. 


APPLICATION OF THEORY TO EXPERIMENT 

Electroepitaxial-growth experiments were carried out 
from Ga-As solutions in a standard LPE apparatus modified 
to permit passage of electric current through the solution- 
substrate interface as described elsewhere.' 5 Semi-insulating 
(Cr doped), nandp type substrates (with a room-tempera-" 
ture free-carrier concentration of about 2x 10 15 /cm : ) were 
used; their thickness ranged from 250 to 1500 pm. 

In the thermally equilibrated growth cell, upon passage 
o current, Joule heating introduced horizontal thermal gra- 
dients (as discussed in detail elsewhere 15 ) and the electroepi- 
taxial growth velocity was found to depend on the solution 
height; this behavior istypical for the presence of convection 
in the solution. The effects of Joule heating were compensat- 
e by modifying the the thermal' configuration of the growth 
system; as the horizontal gradients were drastically reduced 
it was found that the electroepitaxial-growth velocity was 
essentially independent of the solution height, indicating the 
absence of convective interference. The present model will 
e applied to the analysis of the experimental results ob- 


Considering a typical configuration for epitaxial 
growth, as shown in Fig. 1, Peltier cooling at the substrate- 
solution interface is accompanied by Peltier heating at the 
substrate-electrical (back) contact interface. Thus, heat 
transport across the substrate affects the amount of cooling 
at the growth interface; indeed A T p at the growth interface 
. increases with increasing substrate thickness. 3 . 5 

For most practical cases, A T can be approximated as 
follows: 

ATp^Jwpd/lK, ’ ( 20 ) 

where 7 is the current density, v p is the differential Peltier 
coefficient, and at is the thermal conductivity of the substrate 
(see also Appendix B). On the basis of Eq. (20) and taking 
into account the fact that in the solution E/J~p (where p is 
the resistivity of the solution), then Eq. (17a) can be rewrit- 
ten to define a critical substrate thickness d c for which the 

contributions ofthe Peltier effect and electromigration to the 

gre^vth velocity are equal: 


1.0 


g 08 


£ 


o 

.p 

<D 
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I Substrate, n-type, 2x!0 I8 /cm 3 
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FIG. 6 Growth velocity or GaAs from Ga-As solutions at 800 “C as a 
I unci ion of current density; substrate thickness was 300pm. 
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* 1 

tained-both in.the absence and in l he presence ofconvective 
interference 

Absence of convection 

The magnitude of the temperature change in the vicini- 
ty of the interface induced by the Peltier effect was deter- 
mined by positioning a- thermocouple near the substrate-so- 
lution interface (0 3 mm from the interface). It was found 
that fona given cu rrent polarity and at a current density of 10 
A/cm 2 the temperaru re. change ranged from —0.3 to 0.5 °C, 
depending omthe thickness and conductivity type of the sub- 
strate (see Appendix B) Approximate calculations show 
that under the present experimental conditions A T p <AT C , 
where 2lF c is givenin Eq. (19a). Thus, electroepitaxial 
growth in .this case is expected to be dominated by 
electromigration. 

Typical experimental results are shown in Fig. 6,' where 
theelectroepitaxial-growth velocity (as determined by inter- 
face demarcation 23 ) on substrates with a starting thickness of 
300 fim is plotted against current density; the polarity of 
both the n~ and p - type substrate was positive. It is seen that 
the growth velocity on the/?-type substrates (Peltier heating 
at the interface) is only somewhat smaller than on //-type 
substrates (Peltier cooling at the interface). Actually these 
results indicate a Peltier-effect contribution to the growth 
velocity not exceeding about 15%, which is comparable in 
magnitude to the.experimental error associated with the de- 
termination of the growth velocity. 

The linear relationship between the growth velocity and 



FIG 7. Elect roepitaxial-growth velocity of GaAs from Ga-As solutions as a 
function of cur rent density; substrates were //-type 2x 10'Vcm«und-300/im 
'thick. 



FIG. 8 Atom fraction of As in Ga-As solutions calculated from the results 
of Fig 7 and Eq. (6b). Dashed line is the Ga-As liquidus (Ref. 32) 


current density in-Fig. 6 is consistent with the present the- 
ory, e.g., Eq. (6b), since E is proportional to J . Taking 
C X (C S — C,) at 800 °Cc=:4.5 X 10" 2 32 and using Eq. (6b) a value 
of 1.3 X 10' 5 cm/sec is obtained' for the electromigration ve- 
locity of As species, uE, at 10 A/cm 2 ; this value is in agree- 
ment with those reported earlier, 1 . 2 29 

Typical results of experiments carried out at various 
temperatures as a function of current density with/3 5 0-/zm- 
thick //-type substrates (n =2X 10 ls /cm 3 ) are shown in Fig. 
7. It is seen that, consistent with the present model, the. 
gro,wth velocity is linearly dependenton currenPdensity for 
alf temperatures and that, for a given current density, the 
growth velocity increases with temperature. 

From the growth-velocity data of Fig. 7 and Eq. (6b), 
the values of the As concentrations C, (at the growth inter- 
face) were calculated at the various temperatures. These val- 
ues are plotted in Fig. 8, together with those obtained direct- 
ly from the Ga-As diagram 32 (liquidus line). It is seen that the 
experimental results as reduced through Eq. (6b);are in good 
agreement with the Ga-As phase diagram. Thus, within ex- 
perimental error, electroepitaxial growth takes place under 
nearly equilibrium conditions. 

Maximum thickness of epitaxial layer 

Themaximum thickness of electroepitaxial GaAs lay- 
ers was determined for a solution -height A. ranging from 0.3 
to 1 .2 cm at'SOO °C using //-type substrates (// = 2 X 1 0 1 Vcm 3 ) 
300 /zm thick. The time necessary for achieving Z ra iv ranged 
from 1 to 5 h. The results are shown in Fig. 9, where Z niJX //i 
is plotted' against//. The observed linear relationship is .in 
agreement with Eq. (f7a) of the present theoretical 
treatment. 

From the slope d ( Z nux Jr')/dh (Fig. 8) the value of 
juE/D at 10 A/cm 1 is found to be.0.22 cm' 1 . Since 
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FIG. 9. Normalized maximum GaAs layer thickness as a function of solu- 
tion height. Experiments carried out at S00 substrates were n type 
2 X lO’Vcm 1 , 300 fim thick (see text) 


^£=1.3X 10" 5 /sec (as determined above), the value of the 
diffusion constant is found to be D=6X 10 5 cm/sec. This 
value determined on the basis of the present electroepitaxial 
model is in good agreement with the diffusion constant of As 
in an equilibrated Ga-As solution at 800 G reported in the 
literature. 30 

r 

Presence of convection 

Assuming that no thermal gradients are present in the 
growth cell prior to electroepitaxy, they can be introduced 
by passage of electric current causing convection in the solu- 
tion. Thus, horizontal thermal gradients are introduced 
from Joule heating at the asymmetrically located current 
electrodes 15 ; vertical thermal gradients can result from Joule 
heating at the substrate; finally horizontal and vertical ther- 
mal gradients can be introduced at the interlace by tne 
Peltier effect. 

According to the present model of eiecrroepitaxial 
growth, the contributiomof elec Emigration to-the growth 
velocity [second term in Eq. (1 0] is independent of convec- 
tive flow; however, the contribution of the Peltier effect does 
depend on the presence of convection [first term in Eq- (1 1)] 
since the solute boundary layer thickness <5 is a function of 
the convective flow in the solution. 

Electroepitaxial growth of GaAs has beerrcarried out 
under conditions leading to pronounced convection 'em- 
ploying Cr-doped substrates with a thickness of 300 /mi and 
a current density of 25 A/cnr at 900 C. The cui rent-in- 
duced temperature decrease at the growth intei taco A T p w as 
U8*C. Under these conditions the growth velocity was 
found to be dependent on convective flow 1 : the results are 
shown in Fig. 10. For very small values cl solution height 
convection becomes negligible and the growth velocity ap- 


8 (cm) 
Oil 005 



FIG. 10 Growth velocity of GaAs from a Ga-As solution (900 °C, 25 
A/cnr) as a function of solution height. Substrates were Cr doped and 300 
thick v E represents the contribution of electromigration to the growth 
velocity. Estimated values of the solute-boundary-layer thickness <5 are 
indicated. 


proaches the theoretical limit corresponding to the growth 
velocity v E due to electromigration alone. The value of v E 
indicated in Fig. 10 was calculated for the present condition 
from Eq. (10), taking //£- 3.2 x 10' 5 cm/sec for 25 A/cm 3 
(since pE for 10 A/cm 2 was determined to be 1.3 X 10' 5 
cm/sec). Equation (11) can be rewritten as 

§= — AT p'dC\ r \" ' ( 22 ) 

v T C s — C, dTlL 

Thus, knowing v T and AT p the solute-boundary thickness 
can be estimated. Values of <5, estimated from Eq. (22) for the 
results shown in Fig. 10. are indicated in the figure. - 

Thus, the growth behavior presented in Fig. 10 demon- 
strates) in excellent agreement with Eq. (1 1), the enhance- 
ment of thePeltier-effect contribution to the growth velocity 
by convection. 

The contribution of the Peltier effect to the growth ve- 
locity show’n in Fig. 10 (for 300-/1 m-thick substrates) in- 
creases with increasing substrate thickness. According to 
Eqr(21 b), when the thickness of the substrate exceeds a criti- 
cal thickness d c , the contribution of the Peltier effect to the 
growth velocity becomes greater than that of-electromigra- 
tion. Consequently, electroepitaxial growth can take place, 
while electromigration of solute proceeds away from the 
arowth interface, i.e., when a negative substrate polarity and 
a thick p-tvpe substrate (Peltier cooling) are used (opposite 
to that illustrated in Fig. 1). Accordingly, the.behavior en- 
countered in Ref. 21 where growth does take place on thick 
(2 mm) p-type substrates with negative polarity, (opposite to 
that required for solute electromigration towards the sub- 
strate) can now be readily explained. 

SUMMARY AND CONCLUSIONS 

The model of electroepitaxial growth presented here re- 
lates quantitatively the growth behavior to Peltier cooling 
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and electron! igrat ion in terms of fundamental growth pa- 
rameters. The validity of the model was tested with experi- 
mental results on eiectrocnit axial growth of GaAs from Ga- 
As solutions. According to the model, in the absence of con- 
vection in the solution, the growth velocity rs controlled by 
electron) igration and the contribution of the Peltier effect is 
negligible. The dependence of the growth velocity on time, 
current density, and'temperuture, and the dependence of the 
maximum thickness of the epitaxial layer on the solution 
height, given by the model under such conditions, were 
found to be in quantitative agreement with experimental re- 
sults obtained in the absence of appreciable convection. In 
addition, the values of the diffusion coefficient of As in Ga- 
As solutions, of the electromigration velocity, and of the 
equilibrium concentrations of As in Ga-As solutions as a 
function of temperature, calculated from the experimental 
results on the basis of the model, were found to be in good 
agreement with the reported values obtained by other 
means. 

The contribution of the Peltier effect to the growth rate, 
according to the model, can become dominant in the pres- 
ence of convection and is a function of the thickness of the 
solute-boundary layer. In agreement W'ith the.model, the 
growth velocity of GaAs in the presence of convection was 
found to increase with solution height (in the presence of 
thermal gradients convective flow increases and, thus, S de- 
creases with increasing solution height). At very small val- 
ues of solution height the growth velocity approached a val- 
ue which was equal to that calculated from the model taking 
into account electromigration alone. Thus, on the basis of 
the model it is possible to determine the relative contribu- 
tions of the Peltier effect and of electromigration to -the 
growth velocity. It is also possible to estimate the thickness 
of the boundary layer. In addition to the influence of the 
boundary-layer thickness on the contribution of the Peltier 
effect to the growth velocity, the model predicts that the 
Peltier-effect contribution becomes dominant for thick sub- 
strates, consistent with experimental results. 

Electroepitaxy controlled by electromigration (i.e., us- 
ing relatively thin substrates in the, absence of appreciable 
convective flow) takes place under essentially isothermal 
conditions. As expected, under such conditions, epitaxial 
layers exhibit improved surface morphology and decreased 
defeetdensity, and should lead to improved device 
structures. 
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APPENDIX A: SOLUTION OF EG. (2) 

The solution of Eq. (2) for C (x f t ) with boundary condi- 
tions (2a)-(2e) can be obtained as a superposition of solution 
C (£,/ ), corresponding to the case where ail the diffusing so- 
lute is initially contained in a plane at a distance £ from plane 


SOLID 


S0LUT ' 0N P M56 ls 

n 


-|^EI 


x=o 


-I K 
4Z 


AC S C 


FIG 1 L Schematic representation of the change tr* solute distribution b\ 
diffusion and by the presence of an electric field; the solute was originally 
confined in plane 4£(see text) 


x :i (Fig. 11): 



(x-fiEt) 2 \ 
2 (Dt) i/2 J 


where A is an arbitrary constant. 


Accordingly, C (f,/) can be expressed as 


C(£0=”-exp( 


2{Dt) W2 ) 


(Al) 


(A2) 


and; thus, C (x,t) can be obtained as 

C(x,t)-C 0 =A (0 f“c(£0</f+2*<0- (A3) ' 

Jx 

Substituting in Eq. (A3) C (£,t ) from Eq. (A2), one obtains 

C (x,t ) — C 0 =A '{t ) erft / W ftX <A4) 

\2{Dt) wl ) 

where erf is 

erfz =(2/V ,/2 ) f exp(— jr dy), erfez =1— erf 2 . (A5) 
Jo 

From Eq. (A4) and boundary conditions (2d) and (2e), 

A'it) ■ 

C, — Cp 

~ erf [ (6 —fiEt )/2{Dt ) 1/2 ] -erf [ -/nEt/,2(Dt ) 1/2 ]’ 

(A6a) 

B’{t)= ~A\t) erfc [ (<5 —fiEt )/2(Dt ) ,/2 ] , (A6E) 

and thus, Eq. (A4) can be rewritten as 

C(y)-C 0 . 

C.-C. 

_ erf )/?(Pr ) l/2 j -erf [(x-/MEj )/2(Dt ) 1/2 ] 

. erf [(<5 -/xEt )/2{Dt ) 1/2 ] -erf [ -juEt/ 2{Dt ) 1 1/2 ] 

(A7) 

and C (x.t )—C„ forx> S. 

By differentiating Eq. (A7) and taking into account that 
erfez = 1 —erfz, one obtains 

dC _ C 0 -C, 

dx x—o (. Dtrr ) 1/2 
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expf — {uEif/ADi ] _ 

ecfc[ -~jii£t/ 2{Dt) x/1 \ —erfc[(<5— )/l{Dt ) ] ' 2 ] 


APPENDIX B: PELTIER EFFECT AT (GaAs-Ga- 
As) solution interface 


CAS) 


By substituting Eq. (A8) into Eq (A2) and considering that 
erfcoc = 0, the expression for the grow th velocity is obtained: 
for 5= oc (absence of convection) 

^ d-C, (D ) l/2 tzp[-<MEt) 2 /*Dt] 

C s -Ci\ut) erfc [ —fiEt/ 2(Di ) 1/2 ] 


—fiE 



and for finite <5 (presence of convection) 


(A9) 


As pointed out in this paper, knowledge of the tempera- 
ture change A T p at the interface, due to the Peltier effect, is 
essential for assessing the prevailing conditions in electroepi- 
taxia! growth. For this reason, the experimental determina- 
tion of A T p was undertaken for the GaAs system as a func- 
tion of the type of conduction and carrier concentration in 
the substrate, thickness of the substrate, and current density. 

A standard electroepitaxy apparatus was employed for 
the measurements. A thermocouple was positioned in the 
solution about 0,3 mm from the substrate-solution interface. 
Measurements were carried out with p - and tt-tvpe sub- 
strates with carrier concentrations from about 10 u to 
10 19 /cm\ thickness ranging from 250 to 1500//m, and under 
current densities up to 40 A/cm 2 at temperatures of 800 and 
900 °C 


v- 


Cp — Ci 

C s -Q 



exp [ — {jiEt ) 2 /ADt ] 

erfc[ —fj.Et/2(Dt) i/2 ] -ertc[(S-/iEt)/ 2{Dt) ul ) 



(A10) 


Equations (3) and (4) can be obtained from Eqs. (A9) 
and (A10), respectively, taking into account that 


c-c-H 

Co c :~ dr 


AT n 


— (All) 



The temperature changes A T p were measured by apply- 
ing dc. The contribution of Joule heating was determined by 
applying an equivalent amount of ac power. In most in- 
stances, Joule heating of the substrate was negligible in com- 
parison with the Peltier effect, except when the substrates 
were semi-insulating (Cr doped), and their thickness exceed- 
ed 500 jum. 

The results are shown in Figs, 12-14. In Fig. 12 A T p is 
plotted against current density for «-type substrate material 
of varying thickness; A T p is directly proportional to the cur- 



FIG. 12 Changes in temperature at the Ga As-soluticm interface due to the 
Peltier effect as a function of current density at 800 9 C - substrates were n 
type 2 X I0 u /cm*. 


FIG. 13 Changes in temperature at the GaAs-solution interface due to the 
Peltier effect as a function of substrate thickness of 800 *C under 10 A/cm\ 
Substrates were 300 /itn thick. 
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FIG, 14 Changes in temperature at the GaAs-solution interface due to the 
Peltier effect (800 °C, 10 A/cm 2 ) as a function of substrate-earner concen- 
tration Substrates were 300 /zm thick- 


rent density and is negative for positive substrate polarity; 
for a given current density \ AT p \ increases with increasing 
substrate thickness. 

In Fig. 13 the values of AT are given as functions of 
substrate thickness. The polarity of all substrates is positive; 
note that for thep-type substrates (having an opposite Peltier 
coefficient to that of the //-type substrates), A T p is positive. It 
is seen that | A T p | exhibits approximately linear depen- 
dence on. substrate thickness [see Eq. (20)]. Similar depen- 
dence of A T p on current density and thickness as that in Figs. 
12 and 13 has been reported earlier for //-type GaAs materi- 
al”; differences in the absolute values of A T p in the two stud- 
ies are probably due to differences in the the thermocouple 
location and/or differences in heat dissipation from the 
interface. 

The values of AT p as a function of carrier concentration 
are given in Fig. 14 for 300-/zm-thick substrates and a cur- 
rent density of 10 A/cm 2 . It should be noted that the depen- 
dence seen in Fig. 14 is in qualitative agreement with the 
expected behavior of the GaAs Peltier coefficient.” 

Values of A T measured at 900 °C were, within 20%, the 
same as those in Figs. 12-14. Thus, A T p in GaAs electroepi- 


taxy can be estimated fora given carrier concentatration, 
substrate thickness, and current density from Figs. 12—14 
and reasonable extrapolations. 
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MINORITY CARRIER MOBILITY IN p-TYPE GaAs 

W. Walukiewiczf^ J. Lagowskif^ L. Jastrzebski,^ and H. C. Gatos 
Department of Materials Science “and Engineering 
Massachusetts Institute of Technology . 

Cambridge, Massachusetts' 02139 

Abstract 

Theoretical calculations of electron mobility in p-type GaAs 
were carried out taking into consideration the screening effects 
and all major scattering processes. Calculated values of mobility 
are presented ,as a function of carrier concentration, compensation 
ratio and temperature. The basic differences between minority 
carrier mobility in p-type GaAs and electron mobility in n-type 
GaAs are- pointed out. A practical procedure is also presented for 
the evaluation of minority carrier mobilitj T from available electron- 
mobility data. 


'On leave from Institute of Physics, Polish Academy of Sciences, Warsaw 
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Although minority carrier mobility is recognized as a key electronic 

parameter for device application, theoretical treatments of minority carrier 

(l - ) 

mobility -in GaAs and related compounds have thus far been very -limited. ' 

As a result, it has been a common practice to assume that the electron mobility 

in p-type GaAs is equal to the electron mobility in n-type GaAs with the , same 

( 1 - 3 } 

free carrier concentration. Such practice ignores the presence of addi- 

V « (4) 

tional scattering centers (heavy holes) in p-type material. Furthermore, it 
neglects differences in screening effects and electron statistics in n-type 
and p-type material. 

The present mobility calculations utilize the variational method of 
solving the Boltzmann equation -without involving the relaxation time approxi- 
mation. This method has recently been used to evaluate electron mobility in 
n-type GaAs as a function of carrier concentration and’ compensation : ratio. ^ 
Comparison of these calculations with experimental values of electron mobility 
has made possible, the determination of the compensation ratio in n-GaAs. It 
has also been shown ^ ^ that the compensation ratios, obtained from mobility 
measurements are in good agreement with those determined from free carrier 
absorption measurements . In the present paper .the method used in ref . 7 is 
extended to the calculation of the minority carrier mobility in p-type GaAs. 
This extension involves some modifications and additional assumptions.- First 
it is assumed that the minority carrier statistics in the conduction band are 
given by a nondegenerate distribution function described by the crystal lattice 
temperature. In the case of photoexcited electrons this assumption. requires 
that the electrons become thermalized within the conduction band prior to re- 
combining with holes or being trapped at localized centers. For GaAs, in the 
temperature range presently considered, this assumption is justified, since- the 
effective relaxation time is ■ few orders of magnitude smaller than the minority 


carrier lifetime. 
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In addition to the major scattering mechanisms considered in ref. 7 
namely ionized impurity, optical phonon, piezoelectric v scattering and acoustical 
phonon, electron scattering by heavy holes “will be taken into account. Because 
of their large effective mass, heavy holes will be regarded as fixed scattering 
centers. ^ Thus, the overall concentration of ionized scattering centers is 

_ 4. — 

given by + N^, where is the concentration of heavy holes, and 

are the concentrations of ionized acceptors and donors, respectively. Neglec- 
ting the contribution from light holes and assuming that the minority carrier 
concentration is much smaller than p^ one obtains: - N* and thus 



Another modification of the approach of ref. 7 centers about screening 
effects. For a parabolic heavy hole band, the screening length is given 
by the following expression: 


1/^ = 


= 5.80-10 


13 


, An-^ 3 / 2 T 1/2 F 

GvV T F 


-1/2-Olp) 


(1) 


o 

where m^ is the heavy hole effective mass, £ o is the static dielectric constant 

F (T1 ) is the n-th order Fermi-Dirac integral, r] = E^^/k T is the reduced 
n p P * o 

hole Fermi energy* The reduced screening energy, a, which describes screening 
of the electron- ionized impurity interaction is given by 


» 


a = 


2m Irk T 
e D o 


( 2 ) 


where m is the conduction band electron effective mass. The screening energy 
e- 


a = a e /e 

CO q' oo 


/Q ^ 

of the electron-optical phonon interactions has the following form: 

(3) 

where is the high frequency dielectric constant. The screening effects 
arising from the conduction band electrons and the valence band light holes 
are negligible as the concentrations of these carriers are much smaller than 
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the concentration of the heavy holes. 

Compensation is included in the present calculations by introducing the 
compensation ratio, 0^, defined as 0^ = The hole- concentration enters 

into the calculation through the reduced hole Fermi energy, r) 5 the compen- 

p 

sation ratio affects the total concentration of ionized impurities as follows: 


••N . = 


2P, 


i 1-0 


(4) 


The values of GaAs parameters used in ref. 7 have also been adopted in the present 
case. The hole Fermi energy has been determined from the hole concentra- 

tion taking the following effective masses^or the heavy (m^) and light holes 
(m^) : m^ = 0.54 m Q and m^ = 0.089 m Q '. 

The minority carrier mobility y ^ as a function of hole concentration is 
shown in Fig. -1, together with the electron mobility,.' y ^9., results reported' 
in ref. 7. It is seen that the values for y ^ are lower than those for y ^ 
in n- type -material with the same carrier concentration. Furthermore-, y ^ 
decreases more rapidly with carrier concentration, than y ^ » The .calculated 
room temperature mobility ratios y as a function of carrier concentra- 


tion is given in Fig. 2. It is seen that for high concentrations this ratio is 


as low' as 0.6, which clearly . indicates that the previously reported, procedures 


(1-3) 


for determining y^ ^ from values of .y corresponding to' the same free carrier 

0‘S' 

16 “3 

concentration, are not reliable for carrier concentrations exceeding 10 cm . 

-The calculated temperature dependence of y ^ and y ^ is given in 

S 6 

Fig. 3. For the low concentrations the difference between y ^ and y ^ is 
» e e 

more significant at low temperatures. This behavior is apparently due to an 
increased' contribution of the electron scattering by heavy holes relative to 
that by the lattice. For higher carrier concentrations the screening, effects 


become of greater significance and lead to a temperature dependence of y 


(P) 
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entirely different from that of y g v . Actually, for sufficiently low 
temperatures becomes greater than y . 

The above behavior of the minority carrier mobility can be understood if 
one considers the fact that with an increase in the hole concentration, there 
is: (a) an increase of the total concentration of the scattering centers, and 
(b) ,an increase of the heavy hole screening. At room temperature the contribu- 
tion from (a) is more significant than from (b), and the mobility decreases 
with increasing concentration. At lower temperatures the mean electron energy 
in the conduction band becomes much smaller than the screening energy of the 
valence band holes. Accordingly, screening is strongly enhanced and, conse- 
quently, (b) dominates (a). Thus, in contrast to n-type material, in p-type 
GaAs there is a temperature and hole concentration region in which the electron 
mobility increases rapidly with decreasing temperature (Fig. 3) and also with 
increasing carrier concentration. 

The effect of compensation' on the room-temperature electron mobility in 

p-type GaAs has been found to be slightly weaker than in n-type material. The 

computed values of y with the compensation ratio as a parameter are given 
' e 

15 19 -3 

in Fig* 1 for hole concentrations ranging from 10 -10 cm . A comparison 

of the present results on y with those on y^ n ^ (ref* 7 ) shewed that at 

room temperature., y can be .approximated with y as follows: 

e £ 

y (p) (p, ,6 ) = y (n) (n,6 ) - (5) 

e r h’ p r e n 

taking n = and 0^ = (l-K^)/ (3-0.^ where © n and 0^ are the compensation ratios 

in n- and p-type material, respectively. Thus, for example, U iu uncompen- 
sated p-type GaAs with a hole concentration is equal to the electron mobility 
jh n-type GaAs with an electron concentration n = p^ and with compensation ratio 
0 n = 0.33. This approximation is particularly valid for low carrier concentra- 
tion material where screening effects are- negligible. However, even for the 
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higher concentrations considered in this study the error associated with 
‘ this approximation does not exceed 10%. 

From a practical ..point of view it should be noted that the validity of 
the' -above procedure i-s not limited to GaAs. Expression (5) can be applied to 
anyp-type material with GaAS-like energy band structure and- scattering proper-, 
ties. For instance, in the case- of indium phosphide the electron mobility 
values reported for n-type material ^^,11) can ut ^£ ze( j f or . obtaining- the 
electron mobility values in p-type material. 

In summary, it, has been shown that the minority carrier mobility in p- 
type GaAs. is significantly affected by the presence of heavy holes which. act 
as additional scattering centers. The present analysis shows that, only at 
high temperatures and at low carrier concentrations the minority carrier- 
mobility in p— type GaAs equals the electron mobility -in n— type material. 

For higher concentrations the room temperature minority carrier mobility in 
p— type material can be significantly lower than the electron mobility in n— type 
material. The minority carrier mobility was’ found' to exhibit an anomalous 
dependence on temperature and on carrier concentration. A practical procedure 
was presented which permits the estimation of the minority carrier mobility 
using available data on electron mobility in n— type material. 
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FIGURE CAPTIONS 

Room temperature electron mobility as a function of carrier 
concentration. Solid curves represent minority carrier mobility 
in p-type GaAs for different compensation ratios. Dashed curve 
shows the electron'. mobility in n-type GaAs. 

•Ratio of electron mobility in p-type GaAs, to the electron 

mobility in n-type GaAs, as a function of free carrier 

concentration. 

Temperature dependence of the electron mobility in p-type (solid' 
curves) -and n-type (dashed curves) GaAs for two different carrier 

t ]_g _3 

concentrations:' 10 and 10 cm . 
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DERIVATIVE SURFACE PHOTO VOLTAGE SPECTROSCOPY 5 

'A NEW APPROACH TO THE STUDY OF ABSORPTION 
IN SEMICONDUCTORS : GaAs 

^ 7C % 

J. Lagowski, W. Kalukiewlcz, M. M. G. Slusarczuk and H. C. Gatos 
Department of Materials Science and Engineering 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 

ABSTRACT 

Derivative surface photovoltage spectroscopy was achieved with 
wavelength modulation; it was applied to GaAs and' permitted, in a. 
single experiment, the determination of changes in the absorption • 
.coefficient over an energy range of 0.5 to ‘4.3 eV. .Photoionization 
characteristics of deep levels were determined. All known critical , 
point transitions up to - 4.3 eV were clearly' obtained from', second 
derivative spectra. In addition, previously unresolved transitions. . 
were observed at about 2.6 eV. Oscillatory photovoltage in high - - 
magnetic fields was observed, and it was used to identify the transitions 
in the vicinity of the energy gap. .. . . . .. 
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ORIGINAL PAGE 13 
OF POOR QUALITY 

This communication reports on the absorption changes in GaA-s, over a 
wide energy range, involving photoionization - of deep levels, interface states, 
transitions in the vicinity of the- energy gap and higher energy critical 
points. Simultaneous determination o,f these transitions was made, possible 
with a simple experimental approach combining wavelength modulation and surface 
photovoltage spectroscopy with strikingly' enhanced sensitivity through- utiliza- 
tion of the MOS structure characteristics. 

Surface photovoltage spectroscopy . has been employed primarily for surface 
state studies. Its inherently low sensitivity in all standard experimental 
configurations does not permit derivative 'measurements . Wavelength modulation 
of photovoltaic effects associated with p—n, junctions and metal— semiconductor 

barriers has been reported for indirect gap semiconductors. ^ However, p-n 

(20 

junction studies v have been limited strictly to the indirect absorption edge. 

The metal— semiconductor photovoltage spectroscopy has been extended to the 

(3) 

vicinity of direct edges. Extension- of this approach well above or below 

the energy gap is not possible because of the high recombination velocity of' 

minority carriers at the metal interface and the appreciable photo-injection of 

carriers from the metal to the semiconductor. None of the above limitations 

is encountered in derivative surface photovoltage. 

GaAs, n-type, single crystals and epitaxial layers with room temperature 

carrier concentration 1.0 x 10^ J to 3.0 xrlO^cm ^ and with mobility 6,900 to 
2 

3,000 cm /Vsec were employed ,in the present study. The surface photovoltage 
(illumination- induced change of the- surface barrier-) was -.measured with respect 
to a transparent reference gold electrode, evaporated on an oxide layer (800 
to 1000 A' thick) grown anodically^ on the GaAs. With this insulating layer 
a sensitivity of 50 nV was achieved in the surface photovoltage measurements 
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(i.e., 2 to 4 orders of magnitude higher than previously reported) permitting 
the determination of relative changes in the surface barrier as small as about 
lO' 7 . 

The illumination system consisted of a quartz halogen light source and 

a double prism monochromator with a vibrating slit. The illumination intensity 

and amplitude of the wavelength modulation were experimentally selected at 

levels (small) not affecting the measured spectra; typically AAA q was of the 
-3 

order of 10 . Precautions were taken to eliminate the interference of 

spurious signals (arising, for example, from light intensity modulation). The 
measurements were carried out in the temperature range of 95 to 300 K.' 

A typical subbandgap derivative surface photovoltage spectrum obtained 
with a melt-grown GaAs is shown in Fig. 1. The three-peak structure was found 
to be characteristic of all melt-grown GaAs crystals studied, although the ■ 
relative amplitude of the peaks was found to vary with the supplying, source 
of the crystals. In high purity LPE layers, only the intermediate energy peak 
(~ 0.9 eV) was observed, and its amplitude was -about two orders of magnitude 
smaller. Accordingly, the observed transitions are associated with bulk levels 
rather than with interface states. No photoionization transitions related to 
interface states were observed. Thus, .the surface states in GaAs-oxide inter- 
faces previously reported on the basis- of capacitance measurements^ apparently 
are not optically active. 

By extending the surface photovoltage treatment, of surface states ^ ' to 
the transitions involving . bulk levels (in the surface space .charge region) it 
is, shown that the derivative surface photovoltage (d.V)/dhv), is directly pro- 
portional to the derivative of the photo ionization cross section, (dcr)/(dhv)' 
and thus to the derivative of the absorption coefficient. 
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OF POOR QUALITY 

The basic features of the deep level absorption characteristics (Fig. 1) 
coincide with the extensively studied emission characteristics employing photo- 
luminescence. Thus, the sensitivity of the magnitude of the peaks to the 
crystal growth technique and the doping level are identical' in both cases; 
similarly, in both cases the position of the peaks is not sensitive to the 
nature of the donor dopant. Furthermore, the position of the minimum (1.02 eV) 
and that of the zero value of the derivative (1.23 eV) in Fig. 1 are s imi lar 
to the observed energy positions of photoluminescence maxima. Thus, the two 
high energy peaks j correspond to the same deep levels, previously 
observed by luminescence and attributed to complexes involving gallium vacancies 
and donor impurities; the similarity of the absorption and emission energies 
is significant, as it proves directly that there is no measurable Stokes shift 

associated with these levels, contrary to the literature report on Ge-doped 
( 8 ) 

GaAs, In view of the above, the spectrum of Fig. 1 can be analyzed on the 
basis of the quantum- defect model of an impurity-band photoionization cross 

(9) 

section, modified by the introduction of the*Lorent 2 ian broadening F. - 

; Ttie results calculated from this model are shown in Fig. 1. It is seen 

that they are in good agreement with experiment;. furthermore, the present 

treatment enables the separation of the individual contributions of the three 

deep levels and the accurate determination of the binding energies, E . It 

t 

should be pointed out that with the same set of fitting parameters (given in 
caption of Fig. 1) it was possible to obtain good agreement with experiment 
for all samples studied in spite of the fact .that the relative magnitudes 
of the experimental peaks were quite different. It is thus concluded that 
the optically active deep centers in n— type GaAs can be approximated in the 
quantum defect model as acceptor centers, without invoking configurational 
changes . 
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In the region of bandgap transitions second- derivative photovoltage 
-spectra can be directly obtained. A typical spectrum is shown in Fig. 2.- • 

The photo-generated transition in the vicinity of the energy gap, E q , the 
spin-orbit splitting, E q . + A q , and the higher energy structures, and 

+ A^, are clearly resolved. The energy positions and temperature dependence 
of these critical point transitions are in good agreement with those obtained 
by reflectance measurements. ai>12) i s important to point out that the 

-half-width of photovoltage structures E q and E q + A q at 95- K is comparable 

- 613) 

to that obtained at 4.2 K in high-resolution electroreflectance spectra. 

• rr 

' The structure observed in Fig. 2, designated as E 2 , (2.59 eV at 95 K) 

coincides in energy with the transition originally reported in, ref. 11 but 

( 12 ) 

not confirmed in. any later study. . Since this structure, unlike -the others 
in Fig. 2, is characterized by a minimum-maximum,, sequence, it can .be. tenta- 
tively identified with the M 2 or critical point-. 

In order to establish directly the contribution from F 0 -r, transitions 
to the structure in the vicinity of E , derivative surface photovoltage measure 
meats were carried out in high magnetic fields. The results obtained at room 
temperature are shown in Fig. 3. It is seen that (a) the dominant structure 
in the vicinity of E^ is essentially unaffected by magnetic fields; accordingly 

this -structure is not associated with' Y 0 -F , transitions or with free exciton 
* ■ * ' ' 00 

transitions, but rather with residual impurities, and/or defects. Thus, the 
structure at E q cannot be identified with the energy gap. (b) For higher 

photon energies magnetic field leads to an oscillatory magnetophotovoltage 
with the minima and maxima shifting toward higher energies with increasing 
magnetic field. The quantitative account of the involved Landau levels transi- 
tions is complicated in GaAs by the complex nature of its ’Valence band. However 
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from the convergence of the magnetophotovoltage oscillations at zero magnetic 

field the energy gap at 300 X of high quality GaAs is precisely obtained as 

E = 1.440 - 0.002 eV. • 

g- 

In view of the above, it is not surprising that this value is greater 
than those previously obtained from the derivative spectra in the vicinity of 
E o (1.420 to 1.427 eV) Transitions in the vicinity of E q have been exten-' 
sively studied in GaAs, as they are important in understanding laser action. 
However, due to uncertainties in the identification of transitions at room 
temperature, the energy gap values obtained by various methods range from 
1.33 to 1.445. (15) 

The present results (figures 2 and 3) clearly show that derivative 
surface photovoltage is highly sensitive to band— structure optical transitions. 
This sensitivity can be qualitatively explained by the theory of surface photo- 
voltage arising from trapping of photogenerated holes by surface (interface) 
states. The spectral response of surface photovoltage is determined by 

(_d.p) / (1+oLp) where a is the absorption coefficient. Thus, for small values of 

the minority carrier diffusion length, L , (e.g., 10 -5 cm) the surface photo- 
* P 

voltage is essentially proportional to the absorption coefficient; as L p 

increases the sensitivity of photovoltage to absorption coefficient decreases. 

Consistent with this model the amplitude of the high energy peaks in the spectra 

(as those shown in Fig. 2) was found to decrease (with respect to the low energy 

peaks) with increasing L^. However, the high energy (E^+Ap structure could 

be resolved even for high quality GaAs epitaxial layers (L = 6 pm and aL »1); 

P P 

this result can be qualitatively explained considering that for high values of 

(X (high energy region) the holes are generated in the space charge region where 

an electric field (neglected in treatment of ref. 16) as high as 10 3 to 10 4 V/cm 

opposes their diffusion; accordingly, L becomes a much shorter, field-contracted 

P 

diffusion length. 


QUALITY 
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In '-summary , It was 'show that derivative surface photovoltage spectros- 
copy .constitutes a unique means for derivative absorption .spectroscopy studies 
at energies .below and well above the energy gap of GaAs.. This technique 
should be applicable to uni-polar .semiconductors in general. 
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FIGURE CAPTIONS 


Figure 1. First derivative spectrum of subbandgap surface photovoltage of 

15 -3' 

(100) GaAs (melt-grown) -with n = 6 x 10. cm . Thei contribu- 
tion of the individual transitions is shown by dashed lines. 
Parameters used in fitting procedure: = 1.10; 0.905 and 

0.83 eV below conduction band edge; F - 40, 15 and 40 meV-, 
respectively; the Bohr radius of the shallow acceptor 

a = 3.5 x 10 ^cm; the ionization- energy of shallow acceptor 
s 

•E =24.7 meV was adopted from ref. 10. 
s 


Figure 2. Second derivative surface' photovoltage spectra or (.iuu; GaAs 

I* 

15 -3 -4 

’ (n = 6 x 10 cm , L =* 10 cm, at room temperature) ; dashed 
line 293 K; solid -line 95' K. The amplitude of the peaks for 

' f ' If 

E q + A q , E 2 , E^ and E^-+ been multiplied by a factor 

of 50, 20, 3 and. 3, -respectively. 


Figure 3. Second derivative surface photovoltage of (100) GaAs epitaxial - 

15 —3 

layer (n =10 cm ) taken in magnetic fields at 300 K. The 
direction of the magnetic field and that of the light propagation 
were perpendicular to the surface. 
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OUTDIFFUSION OF RECOMBINATION CENTERS FROM THE SUBSTRATE 
' TO THE EPITAXIAL LAYER IN LPE GROWTH ; GaAs 

* 

L. Jastrzebski, J. Lagowski and H. C. Gatos 
Department of Materials Science and Engineering 
Massachusetts Institute of Technology 
Cambridge,-' Massachusetts 02139 

Abstract 

. Experimental results are presented showing, that outdiffusion of recom- 
bination centers from the substrate to the epitaxial layer takes place during 
growth. Such" out-diffusion decreases the carrier lifetime in the epitaxial 
layer to much lower values than the radiative recombination limit. Furthermore, 
it introduces a lifetime- gradient across the epitaxial layer which depends 
critically on the growth velocity and thermal treatment. High rates of growth 
(such as- those - attainable in electroepitaxy) and high cooling rates- 1 can 
minimize .the adverse -effects of normally available substrates on the epitaxial 
layers; however, good -quality substrates are essential for the consistent growth 
of device quality layers. 


INTRODUCTION 

It has been established that the performance of many semiconductor devices 

is adversely affected by nonradiative recombination centers. Their presence in 

the" device active region leads to a decrease in the device efficiency (as in . 

the' case of solar 'cells ^ ) . Furthermore, the "recombination-ahhanced diffusion 

( 2 ) 

of such centers^ J causes a slow degradation of the device performance (as in 

C3') 

the case of semiconductor lasers v J ) . 


The nonradiative recombination centers in GaAs are usually associated with 

i 

residual impurities such as oxygen, ^ transition metals ^ and nonstoichiometric 
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point defects. The origin, of nonstoichiometric defects is. not' fully under- 
stood, although in GaAs it has been attributed to the solidus of the phase 

diagram^ ^ and/or to As depletion of the substrate during the heating cycle 

(S) 

prior to epitaxial growth. Growth at low temperatures and backmelting 

of the substrate prior to growth v tend to reduce the concentration of non- 
stoichiometric defects and improve the quality of the grown epitaxial layer. 

Poor quality GaAs substrate material is a source of defects which can 
diffuse into the epitaxial layer during the growth process. The available melt- 
grown GaAs substrates are in most instances highly compensated, exhibit 
large carrier concentration inhomogeneities and contain high concentration 
of nonradiative recombination centers. (12,13) 

The present work is concerned with the outdiffusion of recombination 
centers from GaAs substrates to GaAs epitaxial layers. 'Minority carrier dif- 
fusion length and lifetime profiling obtained with EBIC-mode (SEM) measurements 
are employed for the study of such outdiffusion. Growth conditions required to 
minimize the outdiffusion of recombination centers are pointed out. 

' EXPERIMENTAL 

Epitaxial Growth and Sample Preparation 

LPE growth of GaAs was carried out employing classical thermal cooling 

techniques and electromigration-controlled electroepitaxy. 14,15 ^ Cd-doped 
17 -3 

(p -' 2 x 10 cm ) 300 ym thick (100) substrates were used. Epitaxial growth 
was performed in a two-well graphite boat. In one of the wells the substrates 
were backmelted (about 25 ym) with under-saturated solution to remove the As 
depleted layer formed during the heating cycle. After backmelting the substrate 
was positioned in the second well containing 2.5 g of undoped Ga— As solution 
with a GaAs source at the top. In the electro epitaxy experiments layers ranging 
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9 

in. thickness from 10 to 100 ym were grown .at 900°C on an area -of 0.5 cm“ 

2 

at a rate of about 6 ym/min by passing an electric current of 60 A/cm for 
a period of 2 to 20 min. Growth was terminated by turning the current off. 

In the thermal growth experiments the growth was performed from equilibrated 
solutions by cooling from 910° C to 900°C with a rate of ,l°C/min. 

In 'both procedures, after growth was completed, the system was kept at 
the growth temperature for time periods ranging from 1 minute to 3 hours with 
the grown 'layer in contact with the solution. The -experiment was terminated 
by quenching the system to room .temperature with an initial cooling rape (from 
900 c C to 700°C) of about 70°C/min. The solution was left over the epitaxial 
layer or wiped by the moving slider prior to cooling. 

„ After, cooling,, if the excess Ga-GaAs mixture was not wiped, it was removed 

from the epitaxial layer -surface by-* boiling in*9Cl; th&- epilayer- was cleaved 

* 

and etched in AB etchant for about 30 sec. Interference contrast microscopy 
was employed to determine the thickness of the epitaxial layers. Whenever 
necessary, the part of the epitaxial layer grown during cooling to room temper, 
ture (abo.ut 5-10 ym in thickness) was etched away in HjSO^ftjO: :H202. Ohmic 
contacts were soldered on the substrate and the epilayer employing Sn and In 
in an H 2 atmosphere. 

. . . The ^epitjjxial -layers -were n- type.- with a carrier concentration of about 
16 -3 

5* X- 10 cm . Schottky barriers required for EBIC-mode measurements were 
made by evaporating aluminum or gold on the surfaces of the epitaxial layer. 
Minority Carrier Diffusion Length and Lifetime - 
. Measurements of .the minority carrier diffusion length were performed 
utilizing electron-beam excitation as shown schematically in Fig. 1. The- 
.diffusion length was obtained from the dependence’ of the electron' beam:- induced 
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current (EBIC) on the distance between the generation position and the 
collecting Schottky barrier or the p-n junction. A 35 keV electron beam 
energy was used in all experiments to minimize the effect of- surface recomr 
bination on the measured .diffusion length. 

The configuration shown in Fig. 1 made possible the determination of the 
diffusion length, L , of -the minority carriers in the epitaxial layer near -the 

Jr 

surface of the epilayer surface (with the electron beam at position 1 and 
with current circuit I) , and near the substrate (with the electron beam in 
position 2 and with current circuit I) ; the minority carrier diffusion length 
was also determined in the p-type substrates (with the electron beam in position 
3 and with current circuit II) . 

Values .of lifetime, x, were estimated from the standard relationship, 

* L - (~ Ut) , where y is the minority carrier mobility and kT/e equals 0.026V 
at 300 K. In p-type GaAs (substrates in- the present case) the mobility of 
electron can be noticeably smaller than the mobility of electrons 
in n-type material of similar free carrier concentration. Accordingly, 

.in evaluating the electron lifetime in the substrates recently calculated theoreti- 
cal values of electron mobility in p— type GaAs were used. In these calcu- 

lations the contribution of electron scattering by heavy holes and the difference 
in screening energies between holes and electrons have been taken into account. 

In the case of holes their mobility as majority or minority carriers is expected 

to be similar at room temperature and for free carrier concentrations below 
17 -3 

10 cm (since screening effects, and carrier- carrier scattering are not signi- 
ficant). Accordingly, in evaluating the hole lifetime in epitaxial layers, 
the hole mobility value of y^ - 300 cm /V-sec was used, i.e., the- hole mobility 
in p-type GaAs with hole concentration of 5 x lO^cm "V^^) which is similar 
to the electron concentration in the epitaxial layers. 
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RESULTS AND DISCUSSION 
Substrate ■ 

The' minority- carrier lifetime in the substrate material was found to be 
about 10 ^ sec, i.e., two orders of magnitude smaller than the lifetime value 
expected from, band-to-band radiative recombination. it was also determined 

that the substrate material' is. highly Inhomogeneous with local lifetime- fluctu- 
ations exceeding a factor of 4. These findings clearly indicate that the life- 
time in the substrates is controlled- entirely by nonradiative -recombination and 
thus' the substrate must contain a high, nonuniform, concentration of- recombina- 
tion centers. 

‘ Outdiffusion of Recombination Centers from the Substrate 

The lifetime in the epitaxial layers was measured as a function of posi- 
tion on a given plane parallel to the surface, as a function of distance from 
the- - substrate-epitaxial layer interface -and as a function of time the layer 
was kept at the growth temperature. 

Typical results of the lifetime measurements as a function of position on 
a plane parallel to the surface of the epitaxial layer are given in Fig. 2. 

The lower curve - corresponds to positions near the surface of the epitaxial 
layer. The layer was 40- pm thick; it was grown electroepitaxially and was 
kept 20 minutes at the growth temperature (from the beginning of growth to the 
beginning of cooling). 'It is' evident from Fig. '2 that lifetime inhomogeneities 
of a similar nature are present near the substrate-epitaxial interface and near . • 
the surface of the layer; however, near the .surface the -lifetime values are 
significantly higher. - It should be noted that -even the highest values of life- 
time are- well below the. values (of the order of 10 ^s) expected from radiative 
band-to-band recombination. 
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The above behavior of the minority carrier lifetime can be understood 
if one considers that the substrate represents a source of recombination 
centers which diffuse into the epitaxial layer during the .growth process and 
thermal treatment* On this basis, the difference between the lifetime near 
the epitaxial layer surface and the lifetime near -the epitaxial layer*- substrate 
interface should decrease with increasing exposure time of the -system to the 
growth temperature. Similarly, the lifetime in the epitaxial layer should 
decrease with increasing exposure time to the growth temperature* As is seen - 
from Fig, 3, such behavior is indeed found in electroepitaxially and in thermally 
grown layers. The lifetime in both electroepitaxially and thermally grown 
layers behaves similarly upon exposure to high temperature after growth, as ‘ 
it is dominated by recombination at outdif fused centers. Without exposure to 
high, temperature, after growth, the lifetime in electroep it axial layers is 
significantly higher than in those grown thermally, as will be discussed below. 

The fact that recombination centers outdiffuse from the substrate into 
the epitaxial layers is further supported by the results shown in Fig. 4. Here 
the lifetime measured near the epitaxial layer surface is plotted as function 
of thickness of the epitaxial layers, all of which have been exposed to a high 
temperature (growth temperature) for the same period of time (20 minutes). The 
layers were grown by electroepitaxy and at the same growth rate of about 6 jim/min. 
It is seen that the lifetime increases by about* two orders of magnitude (near 
the surface) as the distance from the substrate (thickness of the layer) increases 
from 5 to 90 pm. 

If the simplified assumption is made that the lifetime is inversely pro- 
portional to the concentration of the recombination centers, the data of Fig. 3 
and 4 can be used to estimate the diffusion constant of the recombination centers. 
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Thus, treating the substrate-epitaxial layer interface as a limited diffusion 
( 21 ) 

source, the lifetime becomes 

1 1 , X 2 . 

t ~m exp (_ 4ot^ a) 

where x is the distance from the substrate-epitaxial layer interface, t is 
the time, and D is the diffusion constant of recombination centers. From this 
expression the solid line in Fig. 3 is obtained by taking x — constant and 
the solid line in Fig. 4 by taking t = constant. From these plots a value for 
the diffusion constant of the recombination centers of approximately 5 -x 10~ 9 cm 2 /se 
is obtained (900 C) . The same value of 5 x 10 ^cm 2 / sec has been reported for 
point defect diffusion in GaAs (gallium vacancies) at 1000°C. ^ 22 ^ 

^ ^- s important to note that the observed lifetime behavior which is con- 
sistent with the above diffuszon model, cannot be explained on the basis of 
impurity segregation effects. A recent experimental and theoretical analysis 
of impurity segregation in GaAs j during electroepitaxy has shown that the maxi- 
mum changes of dopant segregation coefficient do not exceed 40%, i.e., they are 
two orders of magnitude smaller than the changes required to account for presently 
reported lifetime behavior. ^ 22 ^ 

Growth Velocity and the Effects of Outdiffusion 
The diffusion constant of impurities and point defects decreases exponen- 
tially with decreasing temperature, and thus lowering the growth temperature 
should drastically reduce the effects of outdiffusion from the substrate. How- 
ever, there are limitations to the lowering of the growth temperature as the 

attainable growth rate decreases significantly .and single, crystal growth becomes 
problematic. 

On the other hand, it is evident that the effects of outdiffusion can be 
.reduced if the growth velocity, R, is much greater than the velocity of the dif- 
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fusion front propagation v^. In a diffusion process a constant concentration 

dx const 

profile can be approximated as = const.; consequently, v^ = ~ - 


Thus, the condition R » v^ can be expressed as: 


a » (f ) 1/2 


(2a) 


or 


R » ■? 

Q 


(2b) 


Expression ‘2a represents the case where growth is performed for certain time t 
and 2b represents the case where the growth of a layer of a thickness d is 
required. 

In the present case the outdiffusion constant of recombination centers 
—9 2 

was found to be 5 x 10 cm /s. Thus, according to 2b, for a layer 40 pm thick 
the. growth- rate required to reduce significantly the effect of outdiffusion is 
R » 0.75 ym/min. Thus, in the present experiments the growth rate was about 
6 pm/min, which is high enough to satisfy the above condition for reducing the 
significant outdiffusion of recombination centers. Accordingly, it is not sur- 
prising that, without a post growth heat treatment epitaxial layers 40 pn thick 
were consistently grown with a lifetime of about 10 ^ sec. 

The growth velocity in thermal LPE is limited by the slow solute transport 

to the growth interface and by the fact that the relatively small supercooling 

(24) 

must be used to prevent spurious nucleation in the solution. . These limita- 
tions are overcome in electroepitaxy. 

Actually employing electroepitaxy a number of GaAs layers were grown at 
900°C-with rates as high as 20 ym/min. Under such conditions the adverse 
effects of outdiffusion from the substrate were essentially eliminated. 


SUMMARY AMD CONCLUSIONS 

It was found that outdiffusion of recombination centers from a substrate 
to the epitaxial layer takes place during liquid phase epitaxial growth. From 
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the study of the time dependence of the diffusion process a diffusion constant 
—9 2 

of 5 x 10 cm /sec at 900°'C was obtained for the outdiffusion of recombination 
centers. This value is .the same as that reported for the diffusion constant 
of Ga vacancies in GaAs at 1000° C. 

The recombination of minority carriers at the outdif fused defects consti- 
tutes the limiting factor for' the minority carrier- lifetime in- -the. epitaxial 
layers. It has also been shown that the substrate-epitaxial layer interaction 
during the growth process can be reduced by increasing the growth velocity. 

\ I 

Thus utilizing high growth velocities attainable in electroepitaxy, an improve- 
ment of as much as two orders of magnitude in the minority carrier lifetime, 
was observed. In the light of the present results, the improved characteristics 
(with respect to degradation) of GaAs lasers grown under high growth velocities^ 25 
can be readily understood. 

Preliminary experiments performed on substrates obtained -from different 
crystals (and different supplying sources) showed that defects with different 
diffusion constants outdiffuse from- the substrates into the epitaxial layers 
during growth. Although detailed knowledge of the substrate growth conditions 
was not available, it seems that factors .such as arsenic pressure and growth 
instabilities due to convection or temperature fluctuations influence the 
density and type of defect formation in bulk GaAs crystals. Apparently, dif- 
ferences in the type and concentration of recombination centers in" Bridgman 
and Czochralski grown crystals are responsible for the observed differences, 
in quantum efficiency between GaAs diodes fabricated on substrates 
obtained by these two techniques and the poor characteristics of GaAlAs lasers 
fabricated on Czochralski grown substrates as compared with those fabricated 
on Bridgman-grown substrates. 
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Although high growth rates can minimize the decrease of the lifetime 

due to the slowly outdif fusing recombination centers (with D of the order • of 
-9 2 • 

10 cm /sec), they cannot influence the effects of fast-diffusing defects 

—5 2 

•(with D of the order of 10 cm^/sec) since the magnitude of the growth rates 
required are not as yet achievable by any epitaxial technique. These types 
of defects must be eliminated from the substrates used in epitaxial growth. 

* i 
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Although ' high, growth rates can minimize the decrease of the lifetime 

due to the slowly outdiffusing recombination centers (with D of the order of 
—9 2 * 

10 cm /sec) , they cannot influence the effects of fast-diffusing defects 

—5 2 

(with D of the order of 10 cm / sec) since the'; magnitude of the growth rates 
required are not as yet achievable .by any epitaxial technique. These types . 
of defects must, be eliminated -from .the 1 -substrates -used in epitaxial growth. 
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FIGURE CAPTIONS 

Figure 1 Schematic representation of the experimental configuration for 

the determination of the minority carrier diffusion length by the 
EBXC method; see text. 

Figure 2 Minority carrier lifetime profiles along two planes parallel to 
the growth interface; o - near the original 'growth interface 
(position 2 in Fig. 1) ; © - near surface of epitaxial layer 
(position in in Fig. !)■. 

Figure 3 Minority carrier lifetime as a function of time at the growth 
temperature (900°C); A, o near the original growth interface; 

6 ~ near the surface of the epitaxial layer; triangles corres- 
pond to layers grown by electroepitaxy and circles, to thermally 
grown layers; all layers were 40 ym thick. 

Figure 4 Minority carrier lifetime near the surface of the layers as a 
function of the layer thickness. 
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